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Integral and derivative dispersion relations 
for pp and pp amplitudeCl 

A. Alkin^ and E. Martynov"'^ 

^ Bogolyubov Institute for Theoretical Physics, 
Metrologichna 14-b, Kiev, Ukraine, UA-03680 

The methods of integral dispersion relations (IDR) and derivative dispersion relations (DDR) 
are applied for analysis of the data on pp and pp total cross sections and ratios of real to 
imaginary part of elastic forward scattering amplitude. The models of pomeron behaving 
as triple, double and simple pole (with intercept larger than one) in the angular momentum 
plane are considered. Predictions of the models are given for the TOTEM measurements at 
7 and 14 TeV. 



I. INTRODUCTION 

At the LHC energies the TOTEM measurements of proton-proton total cross-section and near 
forward elastic proton scattering are in fact the first "measurement of pure pomeron only" . Many 
estimates made in the various models of elastic scattering amplitude show that /- and a;-reggeon 
contributions are negligible at LHC energies and small momenta transferred. The forward real and 
imaginary parts of a hadron elastic scattering amplitudes are not independent because an ampli- 
tude obeys the integral dispersion relation which is originated from their analyticity. Therefore 
it would be interest and of importance to investigate a goodness of the experimental data fit by 
the IDR method. However to reach this we need the reliable model for imaginary part of scatter- 
ing amplitude (or for total cross section) in whole kinematical region starting from the threshold. 
There are a lot of the models describing the high energy data on at and ratio p of real to imaginary 
part of forward scattering amplitude. But they cannot be directly extrapolated to low energies in 
order to explore the IDR. 

In this paper we suggest and realize the procedure which allow us to use the IDR for analysis of 
any high energy model of the forward pp pp scattering amplitude. Concretely we have considered 
and analyzed three pomeron models and predict the values of cross sections and ratios of real to 
imaginary parts of forward amplitude an LHC energies 7 TeV and 14 TeV. The procedure can be 



*Based on the talk given by E. Martynov at the Workshop "Forward Physics at LHC", Elba, May 27-29, 2010. 
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easily extended for other amplitudes, such as Tr^p and K^p ones. 

II. GENERALITIES AND DEFINITIONS 

Amplitude. S-matrix theory postulates that the amplitude of any hadronic process is an 
analytic function of invariant kinematic variables. For the ab — t- cd processes under interest, i.e. 
p{p)p —7- p{p)p, they are 

s = {Pa+Pbf, t = {pa-pcf, u={pa-pdf, s + 1 + u = ml^ + ml + ml + m\ = Am^. {!) 

Crossing-symmetry. Crossing-symmetry means that processes p + p — > p + p (s-channel) 
and p + p — > p + p (t, n-channels) are described by the limiting values of unique analytic function 
j4(s,t,n) but taken in different regions of the variables s,t and u. As only two of three variables 
s,t,u are independent in what follows we write down often A{s,t) instead of A{s,t,u). 

Structure of singularities. Singularities of pp and pp elastic scattering amplitudes at t = 
are shown in Fig{T] They are: i)the branch points at s > Am^ corresponding to the threshold 
energies of elastic and inelastic processes, 

ii)branch points generated by the thresholds in u-channel at s < 0, iii)nonphysical (for elastic pp 
scattering) branch points generated by n-channel states (at 4m^ < < ). Thus for amplitude 
we have the right-hand and the left-hand cuts. 

? 2 

s-plane 4mp 4(mp+ m,^ 

pp 



pp physical region 



pp physical region 



4mp-4m^ 



FIG. 1: Structure of singularities of pp and pp elastic scattering amplitudes 

The physical amplitudes of pp and pp elastic scattering are determined at the upper side of 
the cut from Am^ to +oo, precisely, App(s,t) = lim A(s + ie,t,u) = A^(s,t,u) at s > 4m^, 
App{u, t) = lim A{s, t,u + ie) = A^{s, t, u) at u > Amp. One can derive from the definition that 



£-S>0 



App{u,t) = lim A{s — ie,t,u) at s + t<0. (2) 



Optical theorem. For pp and pp it is read as 



aZ{s) ^a+ = -L-ImApp{s, 0) = ^-^^ ImA+{s, 0), (3) 

^^- = ^ImApp{s,Q) = -^ImA.{s,Q) (4) 
2mpp '^'^ 2qsyjs 
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where p is a momentum in the laboratory system, q, 




- Ami ^±(s,0) = A%{s.{)). 



Polynomial behaviour. It is weh known that the scattering amphtude may rise at high values 
of \s\ not faster than a power of s, i.e. a such N exists that at |s| — )■ oo and Iq <t <{) 

\A{s,t)\<\sf. (5) 

High-energy bounds for cross-sections. Total hadron cross sections behave at asymptotical 
energies in accordance with the well known Froissart-Martin-Lukaszuk bound 

ct(s) < — :rln^(s/so) at s — )• oo, sq ~ IGeV^. (6) 

The last inequality means that \A{s,&)/s^ \ at |s| — )• oo. 

All these properties of elastic scattering amplitude are of importance at deriving the integral 
dispersion relations. 



III. INTEGRAL DISPERSION RELATION (IDR) 

As analytic function of variable s the amplitude A{s^ t) (in what follows we consider forward 
scattering amplitude, t = 0) must satisfy the dispersion relation which can be derived from the 
Cauchy theorem for analytic function . 

Because of asymptotic behaviour of the hadronic amplitudes (|^(s, 0)/s^| — )■ but \ A{s, 0)/s| 
at |s| — )• oo) it is more convenient to apply Cauchy theorem to the function A{s, 0)/((s— so)(s— si)) 
rather than directly to amplitude. Generally, the points sq and si are arbitrary ones but usually 
they are chosen at sq = si = 2mp. 

Deformingintegration contour C in the Cauchy relation for amplitude (more details can be found 
in the books [ij) and neglecting the integral over circle with infinite radius (because \ A{s, 0)/s^| — )■ 



at \s\ — 7- oo) one can write 



A{s,0)= Aiso,0) + is-so)A'iso,0) + ^^ 



r Ds(s' O) , / r Ds{s',0) , / 
J (s'-so)2(s'-s)"* ^ J (s'-so)2(s'-s)"* 



(7) 



/ (s'-/o)2(Xs)^'^' + / (u'-uoY^iJ-u)^''^' 
4m2 4m2 



= A{so,0) + {s-so)A'{so,0) + ^-^ 
where 

Ds{s,t) = ^[A{s + ie,t,u) - A{s - ie,t,u)], Du{u,t) = ^[A{s,t,u + ie) - A{s,t,u - ie)]. (8) 

After some simple transformations one can obtain the standard form of integral dispersion relations 
written in the laboratory system (s = 2mp{E + mp),u = 2mp{—E + m^), point sq corresponds to 
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^n = 



P±(r± = + + — P 

Inipp p TTp 



+ 



E'^{E'-E) E''^{E' + E) 



p dE' (9) 



where A'{z,Q) = {d/ dz)A{z,Q). 

We would like to note that if the difference of cross sections Ao" 



0"+ — (T_ vanishes with the 



energy, i.e. if odderon does not contribute asyniptotically to the forward pp and pp amplitudes 
(it was confirmed by the COMPETE analysis [4] that there are no indications for any visible 
odderon contribution to (Jtot{s) and p{s)), then it can be proved that the integral dispersion relation 
Eq. ([9]) with two subtractions is reduced to the dispersion relation with one subtraction. Thus, if 



0-4 



cr_ at s — 7- oo the dispersion relations for pp and pp amplitudes in terms of the cross 



sections and parameters p can be written in the following form [2] 



P± cr± 



B E 
- + — P 

p up 



0"=, 



E'{E'-E) E'{E' + E) 



p dE' 



(10) 



where S is a constant to be determined within a some specific model or e.g. from the experimental 
data fit. 

Eq. (|10p does not contain a contribution of the nonphysical cut (0, 4(mp — m^)), however, as 
pointed out in [6] at high energies it is small (less a few percents). Therefore, applying IDR for data 
description at ^/s > 5GeV we suppose that the subtraction constant B treated as free parameter 
efficiently mimics also a small nonphysical contribution. 



IV. DERIVATIVE DISPERSION RELATIONS (DDR). 

The starting point to obtain the derivative dispersion relations is the dispersion integrals with 
two subtractions for an even A^'^^ and odd A^~^ amplitudes 

A^^\s,t) = hA^{s,t)±A.{s,t)] (11) 
1—1 ^ 

It was proved in [3] that at high energy the following formal relations are valid 



ReAM(^E,0) -B^ Et^n (f i?^) ImA(+\E,0)/E ^ Ete^n i?^) '^J^^af . 



(12) 



where B and C are originated from subtraction constants in IDR. At high energy the approximation 
tan{{TT/2)Ed/dE) « {7r/2)Ed/dE can be used. 



In the papers 



a] all corrections to the asymptotic expressions Eqs. (|12p were derived in the 



form of nip/E power series. 
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V. DESCRIPTION OF THE utot AND p DATA AT ^/i > 5 GEV IN IDR AND DDR 

If the imaginary part of scattering amplitude (or total cross section) is known there are two 
ways to calculate its real part or ratio p. 

1) . i?eA(s, 0) is calculated through IDR, however in this case the imaginary part of amplitude 
(or total cross section) must be known in whole kinematical region, just from the threshold. 

2) . ReA{s^^) is calculated through DDR, this method is good for high enough energy, but for 
low energy it is necessary to know many derivatives of scattering amplitude in order to determine 
corrections. 

We compare the both methods fitting the data P] on atot and p at energy y/s > 5GeV in three 
models for pomeron. 

Low energy cross sections in IDR. As it was mentioned above the imaginary part of the 
amplitude under consideration must be known in whole kinematic region for a given process starting 
from the threshold. In our opinion the best way to solve the problem phenomenologically is the 
following. The cross-sections at low energies (153 pp points and 385 pp points at ^/s < 5 GeV) 
are parameterized as much well as possible. The obtained values of the integral can be used with 
the various high energy models of the total cross sections, because they do not affect high energy 
results. 

The parametrization was taken as follows. For pp cross-section: 

Cpi +gpipe~P/Pp\ P <Pi, 

Cp2 + 5p2(l + 9p3eP/P''')/il + gpAeP/P"^), pi< p < P2, (13) 
Cp3 + gpbs"" + gpes^p, p2 < p < Pm- 



For pp cross-section: 

^Zip) = 

In the above expressions p is a momentum in the laboratory system. The value of pm is the 
momentum at Smin- Parameters Cpk, Cak and some ones of g are not free, they are determined from 
the constraints cjtot{pik — 0) = (^totipik + 0). The values of the rest free parameters are determined 
by a data fit at p < pm under the following constraints. The values of the cross sections (j{pm) 
calculated in accordance with the Eqs. (fT3|) and (fH]l must be equal correspondingly to o-^^^{smin) 
and af^^{smin) given by fitting the specified model for pomeron at high energies s > Smin- One can 
be ascertained that such a tuning does not lead to a big variation of high energy parameters. 



^Zip) = 



Cal + gaie P/P"! + ga2e P^P"^, P < P3, 

(14) 

Ca2 + gaSS'''''' + S'a4S"'"% PS < P < Pm ■ 
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As one can see from the Fig. [21 a quality of the data is far to be good, a data spread is wide. 
There are a few points highly deviating from the main bulk of the data, each of them gives > 
50. If the contribution of only these points to the total is not accounted (but points are not 
excluded from the fit) we obtain /N « 1.5 for the parametrization (jl3|14p . so a quality of the 
fit is quite good. 

Having a parametrization of the cross sections at low energy we are able to apply IDR for 

various high energy models. We perform an overall fit in the three steps. 

The first step. The chosen model for high energy cross-sections is fitted to the data on the cross 
sections only (without p data) at s > s„j,;„. 

The second step. The obtained "high-energy" parameters are fixed. The "low-energy" parameters 
from Eqs. p3ll4p are determined by the fit at s < Smin, provided the crPP{smin) are given by the first step. 

The third step. The subtraction constant B+ is determined by the fit at s > Smin with all other high 
energy parameters being fixed (alternatively all high energy parameters are free, the results below are given 
for this case). 

This procedure allow us to calculate the ratios ppp and ppp at all energies above threshold. The 
results are given below for the specified high energy models. 

High energy. Pomeron models. We consider three models leading to the different asymp- 
totic behaviour of the total cross sections. For each model we investigate how the integral and 
derivative dispersion relations work. We start from the explicit parametrization of the total pp and 
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pp cross-section. Then to find the ratios of the real to imaginary part we apply the IDR as it was 
described above and compare the results for ratios calculated as well making use the DDR. 

The models under consideration are parameterized as the contributions of Pomeron, crossing 
even and crossing odd reggeons 

= ^ CP^E) + R+{E) ± R^{E)) , (15) 

where V{E) is a pomeron contribution while R-^-,R- mimic the contributions of all crossing-even 
(/, 02,... ) and crossing-odd terms {uj,p,... ). Actually and i?_ are very close to /- and oj- 
reggeons which are the most important at the intermediate energies. 

R±{E) = 9±z'^^^'^-' (16) 

where zt = \ cos = 2s/{t — Am^)\{= E/rUp at t = 0) is the "natural" Regge variable with 
'dt being the scattering angle in t-channel |l|. 

The well known parametrization with the replacement E/nip — )• —is/sQ, sq = IGeV^ and 
with asymptotic form of the optical theorem, a = {l/s)ImA{s,0) (it is a standard Regge-type 
parametrization which simultaneously gives the both imaginary and real part of amplitude, we call 
it " — is" parametrization) was considered as well for a comparison. 



a. Simple Pomeron model (Modified Donnachie-Landshoff" model 




VsiE) = Zt {go + 5i^r^°^~'} , ap(0) = 1 + e > 1. 
The standard Regge amplitudes ( " —is" rule) are 

{so/s)APP{s, 0) = igo + iffi(-is/so)"^(°)-i + ig+i-is/sor-^^""^'' ± ^-(-^V^o)"-^")-^ (17) 

b. Dipole Pomeron model 

Vd{E) = Zt {go + gi In zt} . 

{so/s)A%{s, 0) = igo + igi H-is/so) + i5+(-is/so)"+(°)-^ ± g- {-is/sor-^'^-\ (18) 

c. Triple Pomeron model 

VriE) = Zt {go + gi In zt + 52 In^ zt} . 

iso/s)A%{s, 0) = igo + igi H-is/so) + ig2 ln\-is/so) + ig+{-is/sor+^''^~^ (19) 
±g^{-is/sor-^''^-'. 
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TABLE I: The values of the pp and pp total cross sections and ratios of the real to imaginary part of forward 
amplitude at the LHC energies obtained in three pomeron models by IDR, DDR and "—is methods. The 
/dof, where "dof" means degrees of freedom (number of data minus number of free model parameters 

Pomeron model 







Simple pole 






Double pole 






Triple pole 




Fit procedure 


IDR 


DDR 


—is 


IDR 


DDR 


—is 


IDR 


DDR 


—is 


CTpp(mb) at 7 TeV 


95.3 


96.3 


91.2 


90.1 


90.1 


90.5 


93.8 


93.5 


91.0 


Ppp at 7 TeV 


0.137 


0.142 


0.109 


0.105 


0.105 


0.105 


0.125 


0.137 


0.107 


CTpp (mb) at 14 TeV 


107.4 


109.2 


100.2 


98.6 


98.6 


99.1 


104.6 


104.2 


100.0 


Ppp at 14 TeV 


0.136 


0.143 


0.103 


0.098 


0.098 


0.098 


0.121 


0.133 


0.101 


xVdof 


1.104 


1.101 


1.137 


1.108 


1.109 


1.146 


1.103 


1.093 


1.137 



VI. RESULTS OF THE FIT 



Omitting details of the fits (postponing them for a forthcoming paper) we concentrate on the 
main results and conclusions. The quality of the fits is presented in the Table, where ~)(^ jdof = 
/{Nexp.points - Nparameters), and illustrated by Figs. El and H 

As one can see from the Table, all models give good descriptions of the data on atot and p. 
Evidently, the fit by IDR method is preferable. The data on a are described with x^/dof ~ 0.91, 
while the data on p are described less well, with a x'^/^p ~ In our opinion it is occurred 
because of the bad quality of the p data. 

We would like to emphasize the fact that the values of p calculated using DDR are deviated 
from those calculated by IDR even at ^/s ^7 — 8 GeV (see Fig. [H]). It means that in order to have 
more correct values of the p at low energies, one must use the IDR rather than explicit analytical 
expressions from the asymptotic DDR. ^ 

Some difference in predictions of COMPETE [4] for atot and p and those obtained there is 
explained mainly by different sets of data. We considered only pp and pp data while in the data 
on irp, Kp and on other hadrons were used. It would be interesting to extend IDR method for all 
available data. 

Concluding we would like to emphasize that, firstly, the IDR method allows to describe the 
data in a wide interval of energy. Secondly, three considered models having a quite good analytical 
properties predict the values of the total cross sections and ratios in not large intervals, at ^/s = 
14 GeV Gpp is in interval 98 - 108 mb. Nevertheless, if TOTEM measure the total cross section 
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IDR 

Simple pole 
Triple pole 
Double pole 



10000 



FIG. 3; Total pp and pp cross sections in the simple, double and triple pomeron models obtained within 
IDR method 
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(GeV) 



4000 6000 8000 10000 12000 14000 16000 
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FIG. 4: The ratios p for pp and pp in 
three pomeron models at all energies. 



FIG. 5: The ratio ppp at LHC energies, pre- 
dictions of three pomeron models 



with accuracy about 1% one of the models can be preferred. 
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Diffractive cross sections and event final states at the LHC 



Konstantin Goulianos 

^ The Rockefeller University 

We discuss a phenomenological model that describes results on diffractive pp and pp cross 
sections and event final states up to the Fermilab Tevatron energy of y/s ~ 1.96 TeV and use 
it to make predictions for Large Hadron Collider (LHC) energies up to -y/s = 14 TeV and 
asymptotically as ^/s — > oo. The model is anchored in a saturation effect observed in single 
diffraction dissociation that explains quantitatively the factorization breaking observed in 
soft and hard pp and pp diffractive processes and in diffractive photoproduction and low 
deep inelastic scattering. 



I. INTRODUCTION 

As we entered a new energy frontier at the Large Hadron Collider (LHC) with data collected 
at ^ = 900 GeV, 2360 GeV, and 7 TeV from Fall 2009 to Spring 2010, it became painfully clear 
that the Monte Carlo (MC) simulations designed to represent the collective knowledge of the field 
on diffractive cross sections and event final states did not meet the challenge presented to them 
in .his new l,|,he. ene.^ environments The most ccn.oniy n.ed even, generates. PVTH.A fl 
and PHOJET [3], were found to disagree not only with the data but also with each other. The 
latter clearly meant that the two simulations could not both be right. Therefore, an update of the 
MCs was urgently needed. Because of the importance of Minimum-Bias (MB) MC simulations in 
estimating trigger rates, backgrounds, and the machine luminosity at the LHC, a "DIFFRACTION" 



a 



that brought experimentalists and theorists 



workshop was organized at CERN on 7 May 2010 
together to exchange ideas with the goal of producing a reliable MC generator for the LHC. This 
paper is based on a talk I presented at that meeting and an expanded version presented at this 
workshop. 

Diffraction dissociation in pp/pp interactions may be defined by the signature of one or more 
"large" and characteristically not exponentially suppressed ^\ rapidity gaps (regions of rapidity 
devoid of particles) in the final state. The rapidity gap is presumed to be due to the exchange 
of a strongly- interacting color singlet quark/gluon combination with the quantum numbers of the 
vacuum, traditionally referred to as "Pomeron" {P). Diffractive processes are classified as single 
diffraction (SD), double Pomeron exchange (DPE), also referred to as central dissociation (CD), 
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FIG. 1: Non-diffractive and diffractive pp interactions. 

and double diffraction (DD). In pp SDp (SDp), the p{p) dissociates while the p{p) remains intact 
escaping the collision with momentum close to that of the original beam momentum and separated 
from the p [p) dissociation products by a forward gap; in DPE both the p and the p escape, resulting 
in two forward gaps] and in DD a central gap is formed while both the p and p dissociate. The 
above basic diffractive processes are listed below, along with two additional two 2-gap processes 
which are combinations of SD and DD and are indicated as SDD: 



TABLE I: Diffractive cross sections, 
acronym basic diffractive processes 
SDp pp-j>p + gap+ [p-J>Xp], 
SDp pp->- [p->- Xp] + gap + p, 
DD pp ^ [p ^ Xp] -\- gap + [p ^ Xp], 
DPE pp^ p + gap + Xc + gap + p, 

2-gap combinations of SD and DD 
SDDp pp ^ p + ga,p + Xc + ga,p+[p ^ Xp], 
SDDp pp^[p^ Xp]gap + Xc + gap + p. 



Here, Xp, Xp and Xc represent clusters of particles in rapidity regions not occupied by the gap(s). 
The 2-gap processes are examples of multi-gap diffraction, a term coined by this author to represent 
events with multiple diffractive rapidity gaps. A special case of DPE is exclusive production, where 
a particle state is centrally produced, as for example a dijet system or a Z boson. 

Below, in Sec. |TI] (strategy) we outline the method we follow to implement an algorithm 
for a MC simulation, in Sec. |IIl] (CROSS sections and final states) we present excerpts from 
previous papers on total and differential diffractive cross sections and final states, and in Sec. IIVI 
we conclude. 
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FIG. 2: Average multiplicity dN/drj (vertical axis) vs. 77 (horizontal axis) for a process with four rapidity 
gaps, Arii{i ==1-4). 

II. STRATEGY 

A phenomenology that is used to make predictions for the LHC and beyond should be based 
on cross sections and final states that incorporate the current knowledge in the field molded into 
a form that can be extrapolated to higher energies. The issues to be addressed is how to take 
into account saturation effects that suppress cross sections, and what formulas to use for event 
final state multiplicity, pseudorapidity, and transverse energy (-Bt) Ql distributions. In addition, 
the structure of an algorithm for implementing this knowledge into a MC simulation should also 
be addressed. The algorithm must be robust against changes in the collision energy, so that it 
may be equally well applied to simulate collisions at fixed target energies as well as at the higher 
energies of pp and pp colliders and in astrophysics. In this section, we outline a strategy that 
addresses these issues. 

The following input is used for cross sections and final states: 

(i) d?a/dS,dt of the diffractive processes listed in Table [J from the renorm model Q]; 

(ii) fTt(s) from SUPERB ALL MODEL Q]; 

(iii) optical theorem — )• Im fei{t = 0) (imaginary part of the forward scattering amplitude); 

(iv) dispersion relations — >• Re fei (t = 0) , using low energy cross sections from global fit ; 

(v) final states: use "nesting" to describe gap processes, where a nest is defined as a region 
of Ar] where there is particle production, in contrast to a gap region where there are no 
particles 



101, 
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Figure [2] shows a schematic r/ topology of an event with four rapidity gaps and three nests of 
final state particles. The cross section for this configuration is presented in Sec. IIII Al 
We propose the following algorithm for generating final states: 

• start with a pp ^ X inelastic collision at ^/s] 

• decide whether the collision is ND or diffractive based on the expected cross sections; if ND, 
use the ND final state expected at ^/s^, if diffractive, select SDp, SDp, DD, or DPE based on 
probabilities scaled to the corresponding cross sections; 

• for each diffractive event, check whether the region of rj where particles are produced, Ary', 
is large enough to accommodate additional diffractive rapidity gaps: if yes, decide whether 
or not the event will have other gaps within this region, again using probabilities scaled to 
the cross sections, and branch off accordingly; 

• continue this process until the region Ar/' is too small to accommodate another diffractive 



It is important to note that in our definition of a ND collision there are no diffractive gaps 
whatsoever in the final state of the event. In this respect, this definition differs from those of 
"inclusive" or "non-single-diffractive" definitions of ND events used in the literature. 



In this section, we discuss briefly the diffraction dissociation and total cross sections using 
information and/or excerpts from Refs. 



gap. 



III. CROSS SECTIONS AND FINAL STATES 



A. 



Diffractive cross sections 



In Ref. the following expression is obtained for the SD cross section [quoting]: 



asd{s,Arj,t) 
dt dAr] 




1 



X 




(1) 



P9ap{Arj,t) 
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where: 

(i) the factor in square brackets represents the cross section due to the wee partons in 
the ry-region of particle production Ar/'; 

(ii) At/ = Ins-Ar]' is the rapidity gap; 

(iii) K is a QCD color factor selecting color-singlet gg or qq exchanges to form the 
rapidity gap; 

(iv) Pgap{A.7],t) is a gap probability factor representing the elastic scattering between 
the dissociated proton (cluster of dissociation particles) and the surviving proton; 

(v) Ngap{s) is the integral of the gap probability distribution over all phase-space in t 
and Ar]; 

(vi) Fp{t) in Pgap{Ar],t) is the proton form factor Fp{t) = e^°^ ... ; and 

(vii) Cgap is a normalization constant, whose value is rendered irrelevant by the renor- 
malization division by Ngap{s). ... 

By a change of variables from Ar] to using Ar]' = InM^ and Ar] = Ins — InM^, 
Eq. ([TJ takes the form: 



(fa{s,M'^,t) _ r o-o 
.16^ 



dM^dt 



^(^) (M2) 



1 + e 



oo 



e6o 



(2) 



where 6 = 6o + lol In [6 is the slope of the diffractive t-distribution] . Integrating 
this expression over and t yields the total single diffractive cross section, 



s — )• oo „ W>v 
Osd ^ 2cjfP exp 



2a' 



(3) 



The remarkable property that the total single diffractive cross section becomes con- 
stant as s — >• cxD is a direct consequence of the coherence condition required for the 
recoil proton to escape the interaction intact. This condition selects one out of sev- 
eral available wee partons to provide a color-shield to the exchange and enable the 
formation of a diffractive rapidity gap. 



Details are presented in Ref. [12!], where this formulation of the cross section is used to derive 
the ratio of the intercept to the slope of the Pomeron trajectory. Good agreement with the ratio 
extracted from measurements is obtained, providing support for the renormalization approach used 
in the phenomenology. 
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A similar expression may be use for DD, DPE, and multigap processes, as discussed in Refs. 10|, 



For example, the differential cross section for the process displayed in Fif. [2] is derived in 



Ref. [li] as [quoting]: 



gap probability 



where the term in square brackets is the pp total cross section at the reduced s-value, 
defined through ln(s'/so) = J^i^v'n {^^^ foi' each gap) is the QCD color factor 
for gap formation, the gap probability is the amplitude squared for elastic scattering 
between two diffractive clusters or between a diffractive cluster and a surviving proton 
with form factor Fp{t), and Ngap is the (re) normalization factor defined as the gap 
probability integrated over all 10 independent variables ti, rji, rj'-, and Ai] = J2i=i ^Vi- 

The renormalization factor Ngap is a function of s only. The color factors are Cg = 
(N^ — 1)~^ and Cg = 1/Nc for gluon and quark color-singlet exchange, respectively. 
Since the reduced energy cross section is properly normalized, the gap probability 
is (re)normalized to unity. The quark to gluon fraction, and thereby the Pomeron 
intercept parameter e may be obtained from the inclusive parton distribution functions 
(PDFs) 



(4) 



13l |. Thus, normalized differential multigap cross sections at t = may be 
fully derived from inclusive PDFs and QCD color factors without any free parameters. 

The exponential dependence of the cross section on Aiji leads to a renormalization 
factor ~ s^*^ independent of the number of gaps in the process. This remarkable 
property of the renormalization model, which was confirmed in two-gap to one-gap 
cross section ratios measured by the CDF Collaboration (see Ref. jisl), suggests that 
multigap diffraction can be used as a tool for exploring the QCD aspects of diffraction 
in an environment free of rapidity gap suppression effects. The LHC with its large 
rapidity coverage provides the ideal arena for such studies. 

B. The total cross section 

In Ref. [si, an analytic expression is obtained for the total cross section using a parton model 
approach and exploiting a saturation effect observed in the SD cross section. The abstract of 



23 



100 







' > 




< 0.05 




■ 


Albrow et al. 




o 
+ 


Armitage et al. 
UA4 


Standard flux 


• 


CDF 




A 


E710 




X 


Cool et al. / 






/ 

/ 


^ Renormalized flux 

" 


ppy 


/kij^ — — 

I 

"knee" at 22 GeV 


T 



10 100 1000 10000 

Vs (GeV) 

FIG. 3: Total pp/pp single-diffraction dissociation cross section data (sum of both p and p dissociation) 
for ^ < 0.05 compared with predictions based on the standard and the renormalized Pomeron flux (from 
Ref. Q). 

Ref. Ql reads [quoting]: 

The single-diffractive and total pp cross sections at the LHC are predicted in a phe- 
nomenological approach that obeys all unitarity constraints. The approach is based on 
the renormalization model of diffraction and a saturated Froissart bound for the total 
cross section yielding at = (t^/so) ■ ln^(s/sir) for s > sp, where the parameters Sq and 
sp are experimentally determined from the y^-dependence of the single-diffractive 
cross section. 

The following strategy is used in Ref. Q| [quoting]: 

• Use the Froissart formula as a saturated cross section rather than as a bound 
above sp: 

at{s> SF) = at{sF) + ^ ■In'' 

• This formula should be valid above the knee in vs. ^/s at \/^p = 22 GeV 
(Fig. ED and therefore vahd at ^/s = 1800 GeV. 

• Use m? = So in the Froissart formula multiplied by 1/0.389 to convert it to mb~^. 
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• Note that contributions from Reggeon exchanges at ^/s = 1800 GeV are neghgi- 
ble, as can be verified from the global fit of Ref. Q]. 

• Obtain the total cross section at the LHC: 

_ / „LHC „CDF\ 
^LHC ^ ^CDF + :L . ( £ £ ) 

So \ SF SF J 

For a numerical evaluation of a^^^ we use as input the CDF cross section at ^/s = 
1800 GeV, cjf = 80.03 ± 2.24 mb, the Froissart saturation energy = 22 GeV, 
and the parameter Sq- 

...The resulting prediction for the total cross section at the LHC at ^/s = 14 TeV is: 
af/^^v = (80 ± 3) + (29 ± 12) = 109 ± 12 mb. 

For ^/s = 7 TeV, the predicted cross section is: 

alf^% = 98 ± 8 mb [at = 7 TeV] , 

The result for ^/s = 14 TeV is in good agreement with a^"^^ = 114 it 5 mb obtained by the global 
fit of Ref. Q], where the uncertainty was estimated from 5e and the s*" dependence from which the 
value of the parameter So was obtained. 



IV. CONCLUSIONS 



We briefly discuss a phenomenological model that describes available results on diffractive pp 
and pp cross sections and event final states up to the Fermilab Tevatron energy of ^/s = 1.96 TeV 
and refer the reader to previous publications for further details. We also outline a procedure to 
be used to implement the predictions of the model into a Monte Carlo simulation that is robust 
against changes in the collision energy, so that it may be equally well applied to simulate collisions 
at fixed target energies as well as at the higher energies of the Tevatron, the LHC, and beyond. 
The model is anchored in a saturation effect observed in single diffraction dissociation that explains 
quantitatively the factorization breaking observed in soft and hard pp and pp diffractive processes 
and in diffractive photoproduction and low deep inelastic scattering. 
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Problems of phenomenological description of elastic pp scattering at the LHC; 

predictions of contemporary models 



Vojtech Kundrat*, Milos Lokajicek'''^ and Jan Kaspar-f, Jiff Prochazka^ 

^Institute of Physics, AS CR, v.v.i., 182 21 Praha 8, Czech Republic 
^CERN, Geneva, Switzerland 

The standard description of common influence of botii the Coulomb and hadronic elastic 
scattering in the proton - proton elastic collisions at high energies with the help of West and 
Yennie complete amplitude is shown to be theoretically inconsistent. The approach being 
based on the eikonal model amplitude removes these troubles. The preference of its applica- 
tion to the analysis of experimental data and in obtaining the predictions of contemporary 
models for proton - proton high energy elastic hadronic scattering are discussed. 



I. INTRODUCTION 



Differential cross section of elastic scattering of charged nucleons at high energies can be 

defined as j_ _ 

^ = ^\F^^''{s,t)\'. (1) 
at sp'^ 

Here F'^^^ {s, t) represents the complete elastic scattering amplitude which has been decom- 
posed according to Bethe [l| into the sum of Coulomb component F'-^ [s, t) known from QED 
and the elastic hadronic component F^{s,t) bound mutually by a relative phase a\l/(s,i): 

F^+^(s, t) = e^"*(^^*)F^(s, t) + t); (2) 

s is the square of the energy in the center-of-momentum system, p is the momentum of 
incident nucleon in this system and a — 1/137.036 is the fine structure constant. 

The complete (simplified) elastic scattering amplitude F'-'''^^ {s,t) has been proposed by 
West and Yennie 2] (for details see, e.g., jS]) as 

F^+^is,t) ^ ^A(t)/2(t)e-*(^'*) + ^^pV^(p(.) + ^)e^(«)*/2. (3) 

Formula ([3]) is valid provided the hadronic elastic amplitude (the second term on its right 
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hand side) has a constant difFractive slope B together with constant quantity p (the ratio of the real 
to imaginary parts of hadronic amplitude) in the whole kinematically allowed region of t. Together 
with the total cross section atot they may depend only on the energy. The two dipole form factors 
fi{t), f2{t) describe the electromagnetic structure of nucleons. 

The t dependence of the relative phase a'^{s,t) in Eq. ([3]) for pp scattering has been estimated 
by West and Yennie |2] to be 

a^{s,t) = - a{\n.{-B{s)t/2) + -f) (4) 

where 7 = 0.577215 is the Euler constant. 

In earlier analyses the data sets have been usually divided into two regions: \t\ < 0.01 GeV^ and 
higher \t\ ^\. It has been assumed that the Coulomb and hadronic interactions interfered in the 
first region while only the hadronic interactions described with the help of phenomenological elastic 
hadronic amplitude F^{s,t) (having usually more complicated t dependence than the hadronic 
amplitude in Eq. Thus two very diverse formulas for the description of elastic differential 

cross section in the two different regions of t have been used. 

It might seem to be more suitable to use instead of Eq. ^ the integral formula 



^WY{s,t) = - In — + / 



t' -t\ 



(5) 



derived earlier by West and Yennie [4]. However, it may be hardly possible if the rel- 
ative phase a'^wY{s,t) should be real. In such a case the phase C'^{s,t) (defined by 
F^{s,t) = i\F^ {s,t)\e~''''^^^'^^) - or p{s,t) - must be t independent, which is in disagreement with 
experimental data. 

Therefore, the decisive preference should be given to a more suitable approach based on the 
eikonal model formulated within an impact parameter representation of elastic amplitudes. 

II. IMPACT PARAMETER REPRESENTATION OF SCATTERING AMPLITUDES 
AND MEAN VALUES OF IMPACT PARAMETER 



It has been shown in papers of Adachi et al. 



7| that the scattering amplitude -F^(s,t) may be 



related to the eikonal 5(s,6) by the Fourier-Bessel (FB) transformation 

s 

where is the two-dimensional Euclidean space of the impact parameter b. 

If formula ([6]) is to be applied at finite energies some problems appear as the amplitude F^ {s, t) 
is defined in finite region of t only. However, mathematically consistent definition of FB transfor- 
mation requires also the existence of inverse transformation. Thus the amplitude should be defined 
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also in its unphysical region [tJ. The uniqueness of that problem has been established by Islam 
S] by continuing analytically the elastic hadronic amplitude F^(s,t) from the physical to the 
unphysical regions. 

Then the elastic hadronic amplitude in the impact parameter space consists of two terms 

tmin 

hei{s,b) = hi{s,b) + h2is,b)= J dtF^{s,t) Joib^/^)+ J dtF^{s,t)Jo{bV^y, (7) 

^min —CO ■ ■ 

a similar relation holds also for the representation of inelastic overlap function ginei{s,b) [lOt] in 
the impact parameter space. Then the unitarity equation in the impact parameter space can be 
approximated as 

lmhi{s,b) = \hi{s,b)\^ + gi{s,b). (8) 
And the total cross section, integrated elastic and inelastic cross sections may be obtained also as 

oo oo oo 

atotis) = 87r bdblmhi{s,b), aei{s) = 8tt bdb \hi{s,b)f , ainei{s) = Stt bdb gi{s,b). (9) 



The functions Im/ii(s, b) and \hi{s, represent then two main impact parameter profiles (for 
total and elastic processes) and describe the intensity of interactions between two colliding particles 
in the dependence on their mutual impact parameter. 

Expressing this amplitude with the help of modulus \F^ {s, t)\ and of phase C^('5, t) it is possible 
to write for the mean-square value of elastic impact parameter - see Ref. [3] 

Jdt\t\ [A\F^is,t)\y Jdt\t\\F^{s,tr(iC''is,t)]' 
{b\s))^, = 4'^^^ + 4* 



/dt |F^(s,t)|2 J dt\F^is,t)\^ 

= {b\s))^^, + {b\s))^^. (10) 

Introducing further the diffractive slope as 



In- 



da 



dt 



d 



F^(.,i)|, (11) 



\F^is,t)\ dt 



one can derive under the validity of optical theorem for the total mean-square value 

Sirf^bdb &2 lmhi{s,b) 



ib\s)),,, =^ —J- = 2B{s,0). (12) 



It is then also possible to write for the inelastic mean-square value 



a 



,2/ N\ CTelis) ,,2/ XV , CTinelis) ,,2, 



crtot[s) (Ttotys) 
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III. EIKONAL MODEL AND INFLUENCE OF COULOMB SCATTERING 

The complete elastic eikonal equals the sum of the Coulomb and the hadronic eikonals and the 
Eq. Q can be rewritten as 

F^+^(s,t) = F^(s,t) + F^(s,t) + ^ / dVi^^(s,9")i^'^(s,[g-gf). (14) 

ITS J 

Eq. ()14p containing the convolution integral between Coulomb and hadronic amplitudes differs 



significantly from Bethe's Eq. In its final form (valid at any s and t) it equals 



12 



F^+^{s,t) = + —h{t)h{t) + F''{s,t)[l - iaG{s,t)], (15) 



where 





= / {in (9 ^ lMt')Mt')] + ^ [™ - l] /(MO) . (16) 

_4p2 

27r 

I{t,t') = I d^" ^^^^"^^'^^^"\ t" = t + t' + 2Vt¥cos^". (17) 



Instead of the t independent quantities B and p, it is now necessary to consider corresponding t 
dependent quantities B{s,t) defined by Eq. ([TT]). and p{s,t) and (Ttot[s) as 

"<••'> = i;;rTO' <'"(o) = ^i-^"(''.* = »)- (18) 

The form factors fi{t) and /2(i) reflect the electromagnetic structure of colliding nucleons; owing to 
the integration over all kinematically allowed region of t in Eq. (|17p their actual t parameterization 
has been taken from Ref. jl3i] : 

«') = E^. ^-i.^ (19) 

k=i 



and 



where the values of the parameters gk and Wk are to be taken from the quoted paper; for futher 
details see Ref. 



3) 

When the Coulomb part in formula (|15p is taken as known the complete amplitude depends 
in principle on hadronic amplitude F^{s,t) only. Thus Eq. ()15p may be used in two complemen- 
tary ways: either one may test the predictions of different models of high-energy elastic hadronic 
scattering given be corresponding hadronic amplitudes F^{s,t), or one may resolve phenomeno- 
logical t dependence of elastic hadronic amplitude F^{s,t) at a given s (and for all measured t 
values), by fitting experimental elastic differential cross section data with the help of Eq. ([1]) and 

Eqs. (HSD-CZD. 
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IV. EIKONAL MODEL AND ASSUMPTIONS OF WEST AND YENNIE 



The given eikonal model a] 



)proach has been origina' 
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energy of 53 GeV in Ref. (data taken from Ref. 



ly applied to the pp elastic scattering at 



14l|). The modulus and the phase of the 



elastic hadronic amplitude F {s,t) have been parameterized as 

\F^{s,t)\ = (ai + a2t)e''i*+^2t2+M^^(c^+C2t)e'^i*+'^2*'+'^^*' , C^(s, t) = arctan -. (20) 

1 — \t/tdijf\ 



The complete elastic scattering amplitude has been t 



len established with the help of Eqs. (fTC|) -PT |) . 



3 



when some quantities have been limited 



In the following we shall demonstrate the fits 
according the crucial conditions imposed in the approach of West and Yennie. The results are shown 
in Fig. 1 where the fits of differential cross section under different assumptions are given. The value 
of corresponding distribution calculated in the whole measured interval of t has been always 
minimized. First, the contribution of the mere Coulomb interaction corresponds to the dashed 
and dotted line. Second, the contribution of the pure hadronic interaction and corresponding to 
p = const, and B = const, is represented by dotted line. When the Coulomb scattering is added to 
this case the graph with a dashed line is obtained. And finally, when the interference between the 
Coulomb and hadronic scattering with p = const, and B(t) (parameterized according to Eqs. (1201)') 
is demonstrated by the full line roughly copying the data. Compared with the case from Ref. 121 ] 
where both the parameters p{t) and B{t) have been fitted, it gives only a little bit higher value of 
X^. It may be mainly the constancy of B{s,t) that must be denoted as fully unacceptable; even if 
the constancy of p must be denoted as unacceptable, too. 

V. MODEL PREDICTIONS FOR pp ELASTIC SCATTERING AT THE LHC 

The eikonal approach has been made use of in deriving some predictions for the region of LHC 
energy. We have made use of the results of four models, proposed by Bourrely, Soffer and Wu 
[isl . Petrov, Predazzi and Prokhudin [lol . Block, Gregores, Halzen and Pancheri (20] and Islam, 



Luddy and Prokhudin 21[, in which the elastic hadronic amplitudes at individual lower energies 
have been mutually correlated. The eikonal approach has been applied to predict the complete 
elastic pp amplitude at 14 TeV, as the West and Yennie approach must be regarded as irregular 
(see Eqs. (l3l) and (jl])). Some of these results have been included already in TOTEM collective 



mm 



papers [16|, 113] and in 



22]. 



It has been possible to determine total cross section crtot{s), diffractive slope B{s,t) and quantity 
p{s, t) with the help of formulas p!T|) and (fTBj) for each model. The integrated elastic hadronic cross 
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FIG. 3: The diffractive slope predictions for pp FIG. 4: The p{t) predictions for pp scattering at 

scattering at 14 TeV according to different mod- 14 TeV according to different models. 

els. 



sections have been determined by integration of modified Eq. ([T]) containing only F^{s,t). The 

values of atot, (^eU and p(0) are given in Table I. It is evident that the predictions of divers 

models differ rather significantly; the total cross section predictions range from 95 mb to 110 mb. A 

higher value of atot has been predicted by COMPETE collaboration [23]: atot = 111-5 it 1.2 ~^2i 
mb having been determined by extrapolation of the fitted lower energy data with the help of 
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5(0) 


P(0) 




[mb] 


[mb] 


[GeV^2] 




JDULli 1 ^i-j Li OiL. 






90 1 Q 


191 

yj. LZ/ L 


Petrov et al. (2P) 


94.97 


23.94 


19.34 


0.097 


Petrov et al. (3P) 


108.22 


29.70 


20.53 


0.111 


Block et al. 


106.74 


30.66 


19.35 


0.114 


Islam et al. 


109.17 


21.99 


31.43 


0.123 



TABLE I: The values of basic parameters predicted by different models for pp elastic scattering at energy 
of 14 TeV 

dispersion relations technique. The predictions of ^ values for higher values of |f| are shown in 
Fig. [21 they differ significantly for different models. The predictions for the t dependence of the 
diffractive slopes B{t) are shown in Fig. [31 Fig. [H displays the t dependence of the quantity p{t) 
that is not constant, either. 

Fig. [5l shows the t dependence of the ratio of interference to hadronic contributions of the ^ 
for all of the given models, i.e., of the quantity 

^(*) - \F^{s,t)\^ • ^^^> 

The graphs show clearly that the influence of the Coulomb scattering may not be fully neglected 
at higher values of either. It is interesting that at least for small \t\ the given characteristics 
are very similar for all four considered models. The impact parameter representation of elastic 
hadronic amplitude allows then to establish different root-mean-square (RMS) values of 

impact parameters that represent in principle the range of hadronic interactions. Their values 
calculated with the help of formulas (|10p -(|13 p for each of the analyzed models are shown in Table 
II. The values of elastic RMS are in all cases lower than the corresponding values of the inelastic 
ones. It means that the elastic pp collisions would be much more central then the inelastic ones 
similarly as in the case of pp scattering at the ISR energies 
denoted as a 'puzzle', see Ref. [25 1. 



24l | for all these models which has been 



VI. LUMINOSITY ESTIMATION AND ELASTIC SCATTERING AT LHC 



Accurate determination of elastic amplitude is very important in the case when the luminosity 
of the collider is to be calibrated on the basis of elastic nucleon scattering. The luminosity C 
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FIG. 5: The t dependence of the ratio of the FIG. 6: The R{t) quantity predictions for pp 
interference to the hadronic contributions to the scattering at 14 TeV for different models from 



^ for pp elastic scattering at 14 TeV according Sec. I VII 
to different models. 



model 


V< bit > 


V< b'i, > 


V< b'L,i > 




[fm] 


[fm] 


[fm] 


Bourrely et al. 


1.249 


0.876 


1.399 


Petrov et al. (2P) 


1.227 


0.875 


1.324 


Petrov et al. (3P) 


1.263 


0.901 


1.375 


Block et al. 


1.223 


0.883 


1.336 


Islam et al. 


1.552 


1.048 


1.659 



TABLE II: The values of root-mean-squares predicted by different models 

relates the experimental elastic differential counting rate ^^{s, t) to the complete elastic amplitude 
F^+^{s,t) (see Eq. (P) and Ref. [4]) by 

Eq. (|22p is to be valid for any admissible value of t. The value C might be in principle calibrated 
by measuring the counting rate in the region of the smallest \t\ where the Coulomb amplitude 
is dominant. However, this region may hardly be reached at the nominal LHC energy due to 
technical limitations. A procedure allowing to avoid these difficulties may be based on Eq. 
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when the elastic counting rate may be, in principle, measured at any t which can be reached, and 
the complete elastic scattering amplitude F'~'~^^ {s,t) may be determined with required accuracy 
at any \t\. However, in this case it will be very important which formula for the complete elastic 
amplitude will be used. We have studied the differences between the West and Yennie 

simplified formula (see Eqs. Q and dJ])) and the eikonal model (Eqs. ()15p -()17 p ). The differences 
may be well visualized by the quantity 

where F^-j^^ {s,t) is the complete elastic eikonal model amplitude, and F^y^(s,t) is the West 
and Yennie one. The quantity R{t) is plotted in Fig. ([6]) for considered models. The maximum 
deviations lie approximately in the center of interference region where the differences between the 
physically consistent eikonal model and the West and Yennie formula may reach almost systematic 
error of 5 % flQ- 
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Total pp Cross Section at TOTEM 



Jiff Prochazka on behalf of the TOTEM Collaboration^ 
^ CERN, Geneva, Switzerland 

One of the main aims of the TOTEM experiment at CERN is to measure the pp total 
cross section at LHC energies, for which current predictions based on diverse approaches and 
models are quite different. Typical estimations give values in the range 90-130 mb at a center- 
of-mass energy of 14 TeV. Only a precise measurement may give the distinction between 
individual models and TOTEM should provide the value with a precision of 1-2%. To 
achieve this goal TOTEM plans to use the so-called "luminosity independent method" . An 
overview of this method will be given. It will be clearly distinguished between experimental 
input, i.e. what is necessary to measure, and model dependent calculations which are also 
needed. 



I. LUMINOSITY INDEPENDENT METHOD 

The TOTEM experiment at CERN [l-4] (TOTal Elastic and diffractive cross section Measure- 
ment) plans to employ the so-called luminosity independent method to measure the pp total cross 
section cJtot- This method is based on the optical theorem: 

atot = ^ImF,V(s,t = 0) (1) 

where F^(s,t) is the elastic hadronic scattering amplitude; p is the momentum of the incident 
proton, t is four-momentum transfer squared and s is the square of the center-of-mass energy. 
Adding the relation for the elastic hadronic differential cross section 

and the relation between luminosity and total rate N^^^^ = + ^mci 

iatot = iVi?+A^inel (3) 

where A'^ is the elastic hadronic rate and N^^^^ is the inelastic rate, we may derive an expression 
for the total cross section which is independent from the luminosity 



IGtt dt 

O'tot 



t=0 



(4) 



inel 
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where the p parameter is defined as 



^ ReFi;^(.,^ = 0) 



One may similarly derive a formula for the luminosity which does not depend on cxtot- 



L 



^2 



incl/ 



IGvr dTvN 



dt 



(6) 



i=0 

The inelastic rate N-^^^^ is a quantity that can be measured "directly" , while the elastic hadronic 
rate A'^ (~di^) can not be measured directly. The reason is that one measures the complete elastic 
rate given in general not only by the hadronic interaction but also by the Coulomb interaction, 
I.e. we measure iVg^+'^ and not Nf^. It is, therefore, necessary to separate the hadronic scattering 
from the Coulomb scattering. This can be done by employing a formula for the complete scattering 
amplitude F (s,t) which takes into account the interference between the Coulomb and the 

Q 

hadronic interaction. One possibility is to use the simplified West and Yennie formula [01 

where 4>{s,t) = =F [in (— ^) + 7] is the relative phase (7 0.5772 is Euler's constant). The upper 
(lower) sign in the relative phase </> in ([7]) corresponds to the scattering of particles with the 
same (opposite) charges. This WY formula has been derived under the assumption that both the 
quantities p and B are t-independent in the whole kinematically allowed region of t. However, the 
assumption of a constant diffractive slope contradicts experimental data of differential cross section 
at higher |t| where diffractive structure has been observed. For other limitations and deficiencies 
of this formula see Ref . Q] • 

Another possibility to make the separation is to use the more general eikonal formula 

FS^^'i^, t) = ^""(^^ t) + ^""(^^ t) [1 T i^Gi^, t)] (8) 



where 

G(s, t) 



mm 

.Mt")f2{t" 



Ht,t')= I d^"^^^^^, (10) 



27r 



t" = t + t' + 2\/tt' COS <I>" and tmin = —s + 4m , see [6|. The upper (lower) sign in ([8]) corresponds 
again to the scattering of particles with the same (opposite) charges. The important property of 
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the eikonal formula is that it properly combines the Coulomb amplitude F^{s,t) for all values of 
t with an arbitrary hadronic amplitude F^{s,t). While the Coulomb amplitude is assumed to be 
well known and is expressed generally as ib^/i(t)/2(t) where a is the fine structure constant and 
/i,2(i) are proton electromagnetic form factors, the hadronic amplitude is unknown; there are only 
several phenomenological models of F^{s,t). For example, in the simplified WY formula ([7]) it is 
assumed that the hadronic amplitude has purely exponential modulus and constant phase in the 
whole kinematically allowed region of t. 

one has to measure first the differential rate 



To determine the elastic hadronic rate 

C+N 



t=0 



^ — down to small values of make the separation of the Coulomb and the hadronic scattering 
to obtain just the hadronic rate and then extrapolate this quantity to optical point t = 
0. To do the separation and extrapolation one needs to know F^{s,t), i.e. a model of elastic 
hadronic scattering needs to be employed. The luminosity independent method is, therefore, a 
model dependent method. 

The last quantity which we need to know to calculate cJtot given by Eq. Q is the p parameter. 
It may be taken as predicted by the COMPETE collaboration [5] (p = 0.132 at y/s = 14 TeV) or it 
may be calculated on the basis of a model of F^{s, t) which is needed also for the extrapolation of 
elasti(^hadronic rate (more consistent approach) . The values of the p also at ^/s = 14 TeV taken 
from [7] as predicted by different models 9l4l5l| are: Islam et al. 0.123, Petrov et al. (2P) 0.0968, 
Petrov et al. (3P) 0.111, Bourrely et al. 0.121 and Block et al. 0.114. If we take, e.g., the value of 
p predicted by the COMPETE collaboration {p = 0.132), the error contribution from (1 + p^) to 
atot at i/i = 14 TeV assuming the fuh COMPETE error band ^ = 0.33 is ±1.2%. 

Hence the quantities to be measured by TOTEM are the complete elastic rate — ^ — and the 
inelastic rate -/Vj^^^^j. Their measurement method will be the subject of the next section. For a study 
of contemporary models of elastic nucleon scattering and their predictions for the LHC see [Tj]. 
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FIG. 1: Elastic differential cross section at y/s ~ 14 TeV as predicted by various models together with the 
i-acceptance ranges for different optics settings. 



The TOTEM apparatus is placed symmetrically on both sides of the interaction point five (IP5) 
of the LHC. Schematic drawings of the "right" arm of the TOTEM detectors are shown in Figs. [2] 
andO 

To detect elastically scattered protons in the very forward direction TOTEM uses a system 
of Roman Pots (RPs) - movable beam-pipe insertions which are equipped with edgeless silicon 
strip detectors designed by TOTEM with the specific objective of reducing the insensitive area 
at the edge facing the beam to only a few tens of microns (~ 50//m). High efficiency up to the 
physical detector border is an essential feature in view of maximising the experiment's acceptance 
for protons scattered elastically (or diffractively) at polar angles down to a few microradians at 
the interaction point. Currently, two RP stations at ±220 m from the interaction point are fully 
equipped with detectors. 

Two tracking telescopes, Tl (Cathode Strip Chambers - CSC) and T2 (Gas Electron Multiplier 
- GEM), centered at ±9 m and ±13.5 m from the interaction point detect inelastically scattered 
charged particles from the interaction point. The Tl coverage in pseudo-rapidity range is 3.1 < 
\ri\ < 4.7, and for T2 the range is 5.3 < |r/| < 6.5. Both telescopes have 2tt (full) azimuthal 
coverage. Simulations show that about 99.5% of all non-diffractive minimum bias events and 84% 
of all diffractive events have charged particles within the acceptance of Tl or T2 and are thus 



II. 



EXPERIMENTAL METHOD 



A. Experimental set-up 
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triggerable with these detectors (all TOTEM detectors are trigger capable). Thus Tl and T2 allow 
to measure the inelastic rate N-^^^^ For more details concerning the TOTEM apparatus see 



Tl: 18-90 mrad 



IPS 




■ 14 m 



friVI-^^ T2 CASTOR 
^^^^^ . (CMS) 



3^ 



FIG. 2: The TOTEM forward trackers Tl and T2 embedded in the CMS forward region. 
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FIG. 3: The LHC beam line on the "right" side of the interaction point IPS and the TOTEM Roman Pots 
at 147 and 220 m (RP147 and RP220). 

B. Measurement of the elastic rate 



To measure the elastic rate 



C+N 



(from which we may separate the elastic hadronic rate 



diVN 



and extrapolate it to t=0, see Section [I]) one needs proton acceptance at small values of \t\. Two 
optics have been proposed by TOTEM for this purpose. One with /3* = 90 m and the ultimate one 
with (3* = 1535 m foreseen at a later stage. The former one can use the standard injection optics 
and thus easier to commission in early LHC operation. The acceptance of RP detectors at 220 m 
at ^/s = 7 TeV for both optics is in Fig. [H This acceptance has been calculated on the basis of a 
simulation for detectors at distance lOcr + 0.5 mm from the beam (a is size of the beam). From 
Fig. |3]we may see that the minimal achievable value of t is |tmin| ~ 0.025 GeV^ for /3* = 90 m and 

l^minl ~ 8 X 

TeV see With such t-acceptance TOTEM will not be able to reach the so-called Coulomb 
region of t values where the elastic scattering is given practically just by the Coulomb interaction. 
The Coulomb-hadronic region (where neither the Coulomb nor the hadronic interaction can be 



10"'^ GeV^ in the case of /3* = 1535 m. of = For the t-acceptance at ^/s = 14 
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El and also Q 



neglected) will be reached instead, see [jI and also [16| for details. Once is known, the elastic 
rate N^^ can be obtained by integration. Because the elastic differential cross section for small 
values of |t| is relatively very high (see Fig. [T]) just a few days of running are needed to measure 
elastic scattering for small values of \t\. 

To discriminate between different models of elastic scattering which give different predictions 
of differential cross section (mainly for higher values of \t\, see Fig. [1]), it is important to measure 
elastic scattering in the widest possible region. The estimated error contribution from elastic rate 
to the total cross section utot is given in Table H] (for details see [3] ) • 




(a) — 90 m; 50% acceptance: 
0.025 GeV^ < |t| < 0.2 GeV^ 



(b) P* = 1535 m; 50% acceptance: 
8 X 10"'' GeV^ < lt| < 0.2 GeV^ 



FIG. 4: The acceptance of RP detectors at 220 m at -^/i — 7 TeV for two different beam optics. 



C. Measurement of the inelastic rate 

The inelastic rate N-^^^^ will be measured by the Tl and T2 telescopes. The main background in 
the cross section measurement comes from beam-gas events which can be largely rejected by pri- 
mary vertex reconstruction. Studies repeated in [l| show that mainly single and double diffractive 
events cause a major loss in the inelastic rate. The undetected single diffractive events are mainly 
those with very low mass below ~ 10 GeV^, see Fig. O since all their particles are produced at 
pseudo-rapidities beyond the acceptance of Tl or T2. The fraction of events which are not detected 
must therefore be estimated by a model. The error contribution to the total cross section coming 
from the total inelastic rates and dominated by the inelastic trigger losses was estimated to ±1% 
for /3* = 90 m and ±0.8% in the case of /3* = 1535 m. 
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80 100 
M [GeV] 

FIG. 5: The acceptance for single diffractive events as a function of the diffractive mass M at -^i = 14 TeV. 







P* optics 




Quantity 


90 m 


1535 m 


dt lt=0 


Extrapolation of elastic cross-section to t = 

(Smearing effect due to beam divergence, statistical errors, uncertainty 
of effective length L^ff. RP alignment, model dependent deviation) 


±4% 


±0.2% 




Total elastic rate 

(strongly correlated with extrapolation) 


±2% 


±0.1% 


A^inel 


Total inelastic rate 

(error dominated by single diffractive losses) 


±1% 


±0.8% 


P 


Error contribution from (1 + p^) 


±1.2% 




(using full COMPETE error band ^ = 33%) 







Total uncertainty in fjtot ±5% ±1 — 2% 
Total uncertainty in L ±7% ±2% 



TABLE I: Estimated error contributions to the total cross section ctot and the luminosity L at = 14 TeV 
from all quantities entering in Eqs. (0]) and 

III. CONCLUSION 



The determination of the total pp cross section based on the luminosity independent method 
requires from theory a model of elastic scattering (i.e., hadronic scattering amplitude F^{s,t)), 
which can be used both for extrapolation of to the optical point t = and also for calculation 
of the p parameter) and a method for the separation of Coulomb and hadronic scattering. The 



elastic rate 



C+N 



and inelastic rate N- , have to be measured. The total uncertainty in dtot at 



cxiivi iiiv^iciouiv. iuu^ ^'inel 

^/s = U TeV is ±5% for /?* = 90 m and ±(1 - 2)% in the case of the final optics /?* = 1535 m. 
The same quantities which determine o"tot determine also the luminosity L given by Eq. Q. The 
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total uncertainty in L is slightly worse because the total rate A''^ + N-^^^^ in Eq. ([6]) is squared. 
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We describe our fct-resummation model for total cross-sections and show its application 
to pp and pp scattering. The model uses mini-jets to drive the rise of the cross-section and 
soft gluon resummation in the infrared region to transform the violent rise of the mini-jet 
cross-section into a logarithmic behaviour in agreement with the Froissart bound. 



I. INTRODUCTION 

The total cross-section is an observable dominated by very large distances: the so far basic 
unsolved problem of QCD. As such, we do not yet have an understanding from first principles of 
its energy behaviour or process dependence, but models abound. Models are based mostly on the 
optical theorem, analyticity, eikonal representation, Glauber theory, Reggeon field theory. One 
most popular and simple model from Donnachie and LandshofF pj describes all known total cross- 
sections with just two terms, one with a decreasing behaviour, from Regge exchanges, and one 
rising, from Pomeron exchange, namely 



<^ah = YabS-"^ + Xab s' (1) 

This model fits reasonably well all known total cross-sections with two universal powers and dif- 
ferent coefficients for different processes. Apart from the fact that there is some evidence from 
photon-photon scattering that the rise is not universal |j2|], a major objection to this model is that 
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it suffers from violating the Froissart bound. The general consensus is that it may give a good 
description in the intermediate energy region, but not at the extremely high energies where the 
Froissart bound is to be valid. 

It has been known for quite some time that the bulk of the total cross-section is due to semi- 
hard QCD processes [3]. Such model input can be used in the eikonal representation, derived from 
an optical description of scattering. Among the eikonal models, there are the so-called mini-jet 
models, which use PDF calculable QCD mini-jets to describe the rise of the total cross-section 
[4-^, or use QCD inspired formulations, among them the papers by Block et al. Models are 
all in need of some amount of data fitting, since none of them is able to obtain normalizations 
from first principles. Other predictions come from Regge based low energy fits implemented by 
analyticity constraints, and limitations from the Froissart bound. 

An example of such strategy is the one followed by the COMPETE collaboration, which is 
reproduced in PDG, and has the following general expression 8|, 

So s s 

where the last two terms reflect the exchange of different Regge trajectories in different processes, 
the first term is sometime referred to as the "Pomeron" and the squared logarithm reflects a 
geometrical picture, in addition to saturating the Froissart bound. Other fits, such as the one 
by Block and HalzenQ], have an additional term linear in Ins. As useful as fits are, one still 
needs models if one wants to learn something fundamental from total cross-section measurements. 
In particular, since the total cross-section reflects the large distance behaviour of the underlying 
theory, we need models which, albeit phenomenologically, probe the infrared region. 

An overview of proton and photon total cross-section data in the medium to high energy range 
[lo| is shown in Fig. [TJ For the sake of comparison, we have plotted photon and proton data, with 
an ad hoc normalization factor. The full line represents the result from a model which we have 
developed through a number of years, and which probes the large distance behaviour through soft 
gluon fc^-resummation with an ansatz for the effective coupling of gluons with quarks in the very 



low momentum region 
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In the following we outline how the model works, using the case of "pp and scattering. Then 
we show how one can apply the model to pion-proton scattering, to predict the vrp cross-section at 
energies reacheable at LHC through neutron detection in the very forward region. 
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FIG. 1: A compilation of proton and photon total cross-sections, as described in the text. 



II. APPLYING THE fct-RESUMMATION MODEL TO pp AND pp SCATTERING 

Our model blends together a QCD calculated input, given by mini-jet cross-sections, and kt- 
resummation. The model is characterized by three different momentum regions, as follows: 

1- Pt ^ Ptmini foi' parton parton collisions, where the perturbative QCD description leading to 
the mini-jet cross-section is applied, with ptmin ~ 

(1 — 2) GeV kept fixed and independent 

of energy; 

2. Aqcd ^ kt < Qmax ^or single soft gluons emitted from initial state quarks before the hard 
parton-parton collision, through the usual asymptotic freedom perturbative coupling as{kt), 



with 



IM Qmax ~ Ptmin In ^/s/pt„ 



3. kt < Aqcd for ultrasoft gluons in a region which is dominated by a singular, but integrable, 
coupling of the gluons with the emitting quarks, aeffih) ~ k^ '^^ as kt — )• 0. 
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Neglecting the real part of the scattering amplitude, the above modeling for the interactions is 
then input to the following expressions: 



-^ABfi „\ ;^AB /I, „\ , aAB /I, „\AB 



1 K rqmax ju n 



(3) 
(4) 
(5) 



where the impact parameter distribution for the collision is obtained from the normalized Fourier 
transform of the resummed expression for soft gluon emission from the initial state in the collision. 
In Fig. [2] we show the energy behaviour of the mini-jet cross-sections for different Parton Density 
Functions (PDFs) as indicated. The lower cut-off for the parton-parton cross-sections is given by 
Ptmin = 1-15 GeV . With such a value and current PDFs, up to y/s 10 GeV the contribution 
from mini-jets to the total cross-section is quite small and the bulk of the cross-section comes from 



n 



AB 



ft{b, s). This quantity can then be parametrized through the convolution of the form factors of 



the colliding particles j^] or through other appropriate methods 

Through these inputs, we describe the rise of the total cross-section . One can follow the various 
steps which lead to our final result through Figs. [21 [3l H] and [5l 
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FIG. 2: The mini-jet cross-section for different 
energies and different densities. 
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FIG. 3: The maximum single soft gluon trans- 
verse momentum in pp / pp scattering. 



The model we just described can be easily extended to pion processes. In Fig. [6] we show 
the pion exchange process which might make the measurement of irp and vrvr total cross-section 
possible at LHC, through the Zero Degree Calorimeters (ZDC) detecting neutrons very near the 
beam direction. The feasibility and interest in such processes has continued for a long time, and 
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has recently been discussed in 13l4l5l]. While expectation at LHC are strongly determined by the 
many lower energy experiments, up to the Tevatron results at y/s = 1800 GeV, the situation for 
pion processes is very different, as direct measurements from target experiments do not extend 
beyond ~ 40 GeV. Thus, predictions for irp and vrvr total cross-sections suffer from lack of data 
in the high energy region, namely the region where the cross-section starts rising because of the 
onset of parton-parton processes. We have plotted in Fig. [7] a preliminary estimate for np total 
cross-section comparing it with predictions from our model for pp and pp scattering {igI ] . In Fig. [7] 



PRS data for vrp have been extracted by Petrov, Ryutin and Sobol [ij] from earlier measurements 
and obtained using the charge exchange mechanism shown in Fig. [6j A recent discussion [l^] 
indicates that our model satisfies factorization for pion and proton total cross-sections. For all the 
processes considered we have also made a comparison with the BSW model [l3]- The comparison 
with the BSW result indicates a good agreement for the pure proton processes but a different high 
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energy behaviour for vrp and a faster increase at low energy. Given the lack of fixed target data, 
this gives further arguments for the advocated need for future measurements in the high energy 
region. More work is in progress and will be published soon. 



Acknowledgments 

G.P. thanks the MIT Center for Theoretical Physics and Brown University Physics Department 
for hospitality. Work partially supported by the Spanish MEG (FPA2006-05294 and FPA2008- 
04158-E/INFN) and by Junta de Andaluci'a (FQM 101). 



[1] A, Donnachie and P. Landshoff, Phy. Lett. B296 (1992) 227. 

[2] R. M. Godbole, A. de Roeck, A. Grau, G. Pancheri, JHEP 0306 (2003) 061; arXiv: hep-ph/0305071. 
[3] L.V. Gribov, E.M. Levin, M.G. Ryskin, Phys.Rept. 100 (1983) 1-150. 
[4] L. Durand and H. Pi, Phys. Rev. D40 (1989)1436. 

[5] A. Grau, G. Pancheri, Y.N. Srivastava, Phys. Rev. D60(1999)114020; e-Print: hep-ph/9905228. 
[6] For a recent review, see P. Lipari and M. Lusignoli, Phys. Rev. D80(2009)074014; e-Print: 
arXiv:0908.0495 [hep-ph] 

[7] M.M. Block, E. Gregores, G. Pancheri and F. Halzen, Phys. Rev D6O(1999)054024; arXiv: hep- 
ph/9809403. 



50 



[8] C. Amsler et al, Phys. Lett. B667 (2008)1. 

[9] M.M. Block and F. Halzen, Phys. Rev. D73 (2006) 054022; arXiv:hep-ph/0510238. 
[10] A. Grau, R.M. Godbole, G. Pancheri, Y. N. Srivastava Eur.Phys.J.C63{2009)69. e-Print: 
arXiv:0812.1065 [hep-ph] 

[11] R. M. Godbole, A. Grau, G. Pancheri, Y. N. Srivastava, Phys. i2euD72(2005)-76001;arXiv: hep- 
ph/0408355. 

[12] A. Corsetti, A. Grau, G. Pancheri, Y. N. Srivastava, Phys. Lett. B382 (1996) 282. 
[13] B.Z. Kopeliovich, I.K. Potashnikova, Ivan Schmidt , J. Soffer, Phys. Rev. D78 (2008)014034. e-Print: 
arXiv:0805.4534 [hep-ph] 

[14] V.A. Petrov, R.A. Ryutin, A.E. Sobol, Eur. Phys. J. C65 (2010) 637. e-Print: arXiv:0906.5309 [hep-ph] 
[15] A.E. Sobol, R.A. Ryutin, V.A. Petrov , M. Murray, Submitted to Eur. Phys. J. C. e-Print: 
arXiv:1005.2984 [hep-ph] 

[16] A. Achilli, R. Hegde, R. M. Godbole, A. Grau, G. Pancheri, Y. N. Srivastava, Phys. Lett. B659(2008) 
137, e-Print: arXiv:0708.3626 [hep-ph]; A. Grau, R. M. Godbole, G. Pancheri, Y. N. Srivastava, Phys. 
Lett.B682 (2009) 55, e-Print: arXiv:0908.1426 [hep-ph] 

[17] A. Grau, G. Pancheri, O. Shekhovtsova, Y. N. Srivastava, Phys. Lett. B693 (2010) 456, e-Print: 
arXiv:1008.4119 [hep-ph]. 

[18] C. Bourrely, C. Soffer and T.T. Wu, Eur. Phys. J C 28 (2003) 97. 



The ATLAS Forward Physics Program 

Sara Diglio" on behalf of the ATLAS Collaboration. 

" LAL, Univ Paris-Sud, CNRS/IN2P3, Orsay, France 
Abstract 

A brief review of the ATLAS forward detector system is presented. The ATLAS forward physics 
program with early data and himinosity determination are introduced and the relevant analysis 
strategies reviewed. A proposed high luminosity upgrade project is also discussed. 

I. THE ATLAS FORWARD DETECTORS 

The aim of forward physics is to study processes in which particles are produced at very small 
polar angle 9 with respect to the beam. In terms of psudorapidity rj = — ln[tan |] it is possible to 
divide the ATLAS detector into a central part and a number of subdetectors to measure forward 
particle production. 

The central detector consists of an inner detector for tracking purposes (|r/| < 2.5), electromagnetic 
calorimeters to measure the energy of electrons and photons {\r]\ < 3.2), hadronic calorimeters for 
measuring the energy of baryons and mesons (|r/| < 4.9), and a muon spectrometer (Ir^l < 2.7). 
The forward detectors are designed to measure forward particle production. These are ALFA, 
LUCID, MBTS and ZDC: they will be presented in the following sections. A layout sketch of these 
detectors is presented in fig. [TJ 

A. ALFA 

The ALFA (Absolute Luminosity For ATLAS) Roman Pots (RP) are located 240 m from 
the interaction point. There will be two RP on each side, the distance between these two stations 
is 4 m. The system provides a pseudorapidity coverage 10.6 < \r]\ < 13.5. The RP system is not 
fixed relative to the beam. At injection the ALFA detectors are in a withdraw position far from 
the beam. When the beam is stabilized the detectors are inserted to the measurement position 
(within 1.5 mm to the beam). 

The ALFA detectors are primarily designed for measurement of the integrated luminosity and 
elastic scattering at small angles. This will be done by measuring the position of elastic and 
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FIG. 1; The layout of the approved ATLAS Forward detectors. 

diffractive protons that have been deflected at very small angles with respect to the beam and 
that pass through arrays of scintillating fiber trackers ( 20 x 64 fibers in each array). These fibers 
provide a spatial resolution of about 30 ^m.. 

In order to realize such kind of measurement ALFA will collect data in special LHC runs at low 
luminosity (10^'^ cm"'^ ■5^"'^), with high /3* optics and reduced beam emittance. 
The infrastructure and part of the mechanics have been installed in the LHC tunnel, while the 
detectors will undergo a final test beam in summer 2010, and be installed along with the final 
mechanics at the first available shutdown. 

Elastic scattering at small angles and luminosity measurement 

The optical theorem relates the forward scattering amplitude to the total cross section of the 
scatterer. It is usually written in the form atot = 47rlm[/e;(t = 0)] where it can be seen that the 
total cross section {atot) is directly proportional to the imaginary part of the nuclear forward elastic 
scattering amplitude extrapolated to zero momentum transfer {fei{t = 0)). 

For small scattering angle, the four-momentum transfer can be determined by —t = {jpi — po)"^ ~ 
{pO)"^ where p states for the proton momentum of the beam [pi and po refer respectively to the 
incoming and outgoing proton) and 6 for the scattering angle at the interaction point. 
The ALFA detector allows to determine the the total cross section and the absolute luminosity 
measuring the elastic scattering down to such small i-values where the cross section becomes 
sensitive to the electromagnetic amplitude via the Coulomb interference term. Since we know very 
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well the electromagnetic amplitude, it will be possible to add a constraint if the Coulomb region 
will be reached. Taking into account the optical theorem and the Coulomb term, the rate of elastic 
scattering at small t-values can be written in a simplified version as : 

— = Ln\fc + fN\ - L7r{-— + — {i + p)e 2) (1) 

where the first term corresponds to the Coulomb and the second to the strong interaction ampli- 
tude, a is the fine structure constant. Fitting the t-distribution allows to determine the absolute 
luminosity (L), the total cross section {(Jtot), the slope (b) and the ratio of the real over the imagi- 
nary part of the forward elastic scattering amplitude {p) without measuring the inelastic rate. It 
is expected that the absolute luminosity will be determined to ~ 2-3 % accuracy. 



B. LUCID 

LUCID (Luminosity measurement using Cerenkov Integrating Detector) [l', 13] is composed of 
two modules located at 17 m from the interaction point on both side of ATLAS. Each module 
is composed of 20 Cherenkov tubes at the end of which photomultiplier tubes (PMT) are placed 
to collect the light. Each Cherenkov tube is made of aluminium, 15 mm in diameter, filled with 
C4F1Q. It covers a pseudorapidity range 5.6 < \rj\ < 5.9 for charged particles with a Cherenkov 
threshold of 10 MeV for electrons and 2.8 GeV for charged pions. 

LUCID is one of the main ATLAS on-line monitors for instantaneous and integrated luminosity 
measurement. The principle of the measurement is based on the fact that the average number of 
interactions in a bunch crossing is proportional to the number of particles detected in LUCID. 
In order to provide the actual luminosity rather than the change in luminosity, LUCID must be 
calibrated using a known absolute luminosity. At the beginning of LHC running the calibration 
procedure is based on LHC parameters using a so called Van Der Meer or beam separation scans. 
The present accuracy from this method is 11 % but this might be improved in the future [4]. Later 
standard candles like Z-boson production can be used. Here the estimated accuracy is in the 5-8 
% range. Ultimately , LUCID will be calibrated using information from ALFA measurements that 
should allow to reach an accuracy of ~ 2-3 % on the absolute luminosity (see sec. [L 



MBTS 



MBTS (Minimum Bias Trigger Scintillators) [l| are segmented scintillator paddles quite close 
to the beam-pipe. The system consists of 32 scintillator paddles, 2 cm thick, organised into 2 disks, 
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one on each side of the interaction point of ATLAS. The system is placed between inner detector 
and end-cap cryostat and provides a pseudorapidity coverage 2.1 < \r]\ < 3.8 . 
The main purpose is to provide a trigger on minimum colhsion activity during the proton-proton 
colHsions at low luminosities. The apparatus is particularly well suited for the measurement of the 
start-up LHC luminosity. Because of heavy radiation, it is expected that the inefficiency of the 
MBTS will increase after some time of higher luminosity operation. For this reason MBTS will be 
active only during the initial running phase where the average number of interactions per bunch 
crossing is expected to be low. 



D. ZDC 

The ZDC (Zero Degree Calorimeter) [4] is placed at 140 m on both sides of interaction point in 
the TAN region (target absorber for neutrals) between the tubes at the point where the single beam 
pipe splits into two. It is a sampling calorimeter composed of four modules: one electromagnetic 
and three hadronic tungsten/quartz calorimeters. The ZDC is able to measure neutral particles at 
pseudorapidity |r/| > 8.3 . 

At the LHC startup phase, in the early data taking period, the electromagnetic module will not be 
installed and its position will be occupied by the LHCf experiment ^. After initial running LHCf 
will be removed and the full ZDC installed. 

The main purpose of the ZDC is to measure the centrality of the collisions in heavy ion runs. During 
proton-proton collision ZDC is used for beam halo, beam gas suppression, luminosity monitor and 
also as an additional minimum bias trigger. It will also be used to tag diffr active processes. When 
the luminosity of 10'^^ cm^'^s^^ will be reached, the ZDC modules will be removed in order to 
minimize the radiation damage. The ZDC will be reinstalled for heavy ion runs. 



II. FORWARD PHYSICS MEASUREMENTS 



Forward region instrumentation at LHC provides a new window to QCD physics. In this section 
a short review of the ATLAS forward physics measurements that have just started with early data 
taking from December 2009 and that are planned to be done in the future is presented. 
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A. Soft Diffraction 



Single (double) diffraction is a low-t process in which a colour singlet (i.e. Pomeron) is exchanged 
between the two protons and one (both) of the protons breaks up into a dissociative system. 
Diffractive events can be tagged by identifying a rapidity gap between the outgoing proton and 
the dissociative system for the single diffraction (SD) or between the two dissociative systems in 
the double diffraction (DD) case. Both single and double diffractive dissociation have large cross 
section of the order of 10 mb. 

There arc two approaches that will be used to measure soft SD at ATLAS. The first one will 
be focused on the invariant mass of the dissociated system Mx and the fractional longitudinal 
momentum loss $, = (where s is the center of mass energy for the proton proton collisions). It 
is clear thet events with low-^ will be contained only in the forward detectors, whereas high-^ events 
will have activity in many areas of the central detector as well. The dissociative system will be 
identified using the inner detector, calorimeters, LUCID and the ZDC. The last two subdetectors 
and MBTS will be used as trigger. 

The second approach implies the use of the ALFA subdetector. The outgoing proton in SD exchange 
can be tagged and measured using special LHC runs with high-^* optics at a luminosity of 10^^ 
cm~^ s~^. ALFA will be able to measure the longitudinal momentum loss ^ directly using ^ = 
1 — 1^ ( where pi and Po are respecively the longitudinal momenta of the incoming and outgoing 
protons). The resolution of the ^ measurement is approximately 8% for ^ = 0.01, falling to ^ 2% 
for ^ = 0.1. 1.2-1.8 million of events are expected in 100 hours at a luminosity of 10^^ cm~^s~^ 



B. Hard Diffraction 

An interesting measurement is to look for hard scattering events with gap on one side of the 
detector. The aim is to study diffractive parton density functions (dPDF), the ratio of single 
diffractive (SD) di-jets to non diffractive (ND) di-jets and the ratio of double pomeron (DPE) to 
single diffractive di-jets. 

A few thousand SD di-jet events in 100 pb~^ with Et > 20 GeV (after trigger prescale in Level 1 
trigger for low transverse energy jets and gap requirement) are expected. 
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C. Central Exclusive Di-jet Production 



Central exclusive di-jet production (CEP) is defined as the process pp ^ p + jj + p where '+' 
states for a large rapidity gap from the outgoing protons. In these events all of the energy lost by 
the protons goes into the production of a hard central system without any other activity. One of 
the possibility is that this system will be constitued from the Higgs boson: this is why recently 
such kind of process has received a great deal of attention. 

A few hundred di-jet CEP events after trigger and analysis cuts with Et > 20 GeV in 20 pb^^ 
of data are expected. The main background comes from double diffractive di-jet production. To 
overcome this background, the idea is to measure the di-jet mass fraction defined as Rjj = 
where Mjj is the invariant mass of the di-jets and M^aio is the mass of all energy deposit in the 
calorimeter. Typically an exclusive event will have Rjj ~ 1 while inclusive/diffractive events will 
have Rjj <^ 1. 

D. Gaps between Jets 

Gaps between jets arises from a 2 — )• 2 scatter via a colour singlet exchange. The typical 
signature is constitued by two high pt jets separated in the detector by a large pseudo-rapidity gap 
At] > 3. This process has been previously measured at HERA [ol and the Tevatron [10] but due 
to the increase of center of mass energy, an improved measurement should be possible at the LHC. 
ATLAS should be able to reach a gap fraction up toArj ~ 9 — )• 9.5. We aim to first rediscover this 
process before studying the dependence on the gap size and jet energies in the kinematic regime of 
the ATLAS detector. Colour singlet exchange also provides a useful early opportunity to study jet 
reconstruction and triggering in the forward calorimeter. Measurements should be possible with 
10 pb^^ of data. 

III. POSSIBLE FUTURE: THE AFP 

The AFP (ATLAS Forward Protons) [7,] is a project of installing silicon and fast timing for- 
ward detectors at 220 m (AFP220) and 420 m (AFP420) from the ATLAS interaction point for 
measurements at high luminosity. AFP420 will approach the beam down to 5 mm while AFP220 
should reach 2-3 mm. Being so close to the beam at luminosities in the 10^'^ — 10^^ cm"'^ 
range requires very radiation hard detectors. Another important requirement is the very good 
position resolution in order to obtain a high mass resolution of the order of 1-2 % for 420+420 
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and 4-5 % for 420+220 configurations. The position and angular resolutions are required to be 
respectively of the order of 10 fj,m and 1 yurad for both AFP220 and AFP420. In order to satisfy 
those requirements, a 3D silicon detector has been chosen. A precision time of flight (ToF) system 
composed of Cherenkov photon detectors with a resolution of 10-20 ps has been chosen to identify 
the primary vertex. This allows to obtain a large reduction in overlap backgrounds. 
The capability to detect both outgoing protons in diffractive and photoproduction processes opens 
up the possibility for a rich QCD, electroweak, Higgs and beyond the Standard Model experimental 
program. A complete summary of the forward physics program at LHC using the proton taggers 



is described in 



8|. 



IV. SUMMARY AND CONCLUSIONS 



ATLAS has a variety of Forward Detectors (MBTS, LUCID, ZDC, ALFA) that will allow to 
exploit several Forward physics topics and to determine the absolute luminosity with high precision. 
We expect to be able to study soft single, double diffraction and gaps between jets when an 
integrated luminosity of the order of 10 pb~^ will be reached. Exclusive di-jet and single diffractive 
di-jet production measurements will be possible increasing the value of integrated luminosity. After 
2010 ATLAS expect to use the ALFA subdetector to study single diffraction and elastic scattering. 
The elastic scattering measurement is interesting in itself but will also be used for a measurement 
of the absolute luminosity of the experiment. Finally the AFP project will allow to extend the 
ATLAS physics program at high luminosities. 
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The TOTEM experiment at the LHC is devoted to a deeper understanding of the pro- 
ton structure with the main goal being a precise measurement of the total proton-proton 
cross-section at LHC energies. In addition, TOTEM will measure elastic scattering over a 
wide |i|-range and perform comprehensive studies of difFractive processes and forward event 
topologies in inelastic proton-proton collisions. This article reviews the status of the experi- 
ment as well as the expected physics, especially in terms of diffraction, from the data taken 
in 2010 and to be taken in 201 1, both in standard LHC and special dedicated TOTEM runs. 



I. INTRODUCTION 

TOTEM is the only LHC experiment exploring the forward region by charge particle mea- 
surement at pseudorapidities r] larger than three and thus TOTEM physics programme, aiming at 
an deeper understanding of the proton (p) structure, is very different compared to the other LHC 
experiments. The main goal of TOTEM is a precise measurement of the total pp cross section atot 
at 1-2 % level using the luminosity-independent method described in detail elsevi^here [l-^- 

Elastic and diffractive protons scattered at very small angles are measured by TOTEM with 
proton detectors embedded in "Roman Pots" (RP), placed along the LHC beam line on both 
sides of interaction point 5 (IP5). Charged particles produced by inelastic interactions in the 
r/-range of 3.1 < |r/| < 6.5 are measured on both sides of IP5 with full azimuthal coverage by 
the inelastic telescopes Tl and T2 of TOTEM embedded in the CMS end-caps. This will allow 
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FIG. 1: The LHC beam line on one side of interaction point 5 with the positions of the TOTEM Roman 
Pots at 147 and 220 m (RP147 and RP220) indicated. 

TOTEM to determine, in addition to atot^ the elastic scattering cross section over a wide |t|-range 
(~10-3 < |t| < 10 GeV2) and study soft diffractive processes and the charged particle event 
topology in the forward region. These studies will be complemented at a later stage by studies of 
a wide range of diffractive and forward physics topics [J], taking advantage of jm unprecedented 
tracking and calorimetry coverage in rj, in common data taking with the CMS [1|| experiment. 

II. EXPERIMENTAL APPARATUS AND STATUS 



To measure elastically and diffractively scattered protons requires the reconstruction of the 
protons tracks by detectors moved as close as ~1 mm from the center of the outgoing beam. 
This is obtained with two RP stations installed, symmetrically measuring protons circulating both 
clockwise ("sector 56") and anticlockwise ("sector 45"), at distances of ~ 147 m (RP147) and ~ 
220 m (RP220) from IP5 (Fig.[T|). Each RP station is composed of two units at a distance of several 
meters allowing a local track reconstruction. Each unit consists of three pots, two approaching the 
beam vertically from the top and the bottom and one horizontally to complete the acceptance for 
diffractively scattered protons. Each pot contains a stack of 10 planes of silicon strip "edgeless" 
detectors with half with their strips oriented at an angle of +45° and half at an angle of —45° with 
respect to the edge facing the beam. These detectors ^, designed by TOTEM with the objective 
of reducing the insensitive area at the edge facing the beam to only a few tens of microns, have a 
spatial resolution of ~20 /im. High efficiency up to the physical detector border is essential in view 
of maximizing the elastic and diffractive proton acceptances. All detectors for the 220 m stations 
were installed in 2009 and succesfully commissioned in 2010, while those for the 147 m stations 
have been tested and will be installed during the 2010/11 winter shutdown. 

To measure charged particles produced in inelastic proton-proton collisions with a very good 
efficiency, TOTEM has the Tl and T2 telescopes embedded in the forward region of CMS (Fig. [2]). 



61 



IPS 



T1: 18-90 mrad 
T2: 3 - 10 mrad 




FIG. 2: The positions of the TOTEM Tl and T2 telescopes in the CMS forward region. The pseudorapidity 
and angular coverage of Tl and T2 are given. 



The Tl telescope consists of Cathode Strip Chambers (CSC) and T2 telescope of Gas Electron 
Multipliers (GEM) Q]. The ry coverage of Tl and T2 is 3.1 < < 4.7 and 5.3 < \r]\ < 6.5, 
respectively. Each Tl telescope arm consists of five planes made up of six trapezoidal formed 
CSCs with a spatial resolution of ~1 mm. Each T2 telescope arm consists of 20 semicircular shaped 
triple-GEM detectors with a spatial resolution of ~100 fim in the radial direction and a inner radius 
that matches the beam-pipe. Ten aligned detectors mounted back-to-back are combined to form 
one T2 half arm on each side of the beam-pipe. The full T2 telescope was installed in 2009 and 
successfully commissioned in 2010. The full Tl telescope has been commissioned with particles 
beams at the CERN SPS H8 test line and will be installed during the 2010/11 winter shutdown. 
TOTEM has been taking data with all installed detectors (both T2 arms and both RP220 stations) 
since the start of LHC operations in 2009. Data taking settings like the timing and thresholds of 
both detectors has been carefully tuned in order to optimize the running conditions. 

The knowledge of the positions of the RPs w.r.t. beam center is a key issue. Therefore a 
collimator-based RP alignment procedure (see Fig. [3]) has been developed in collaboration with the 
LHC operators taking advantage of the sharp beam edges produced by the LHC collimation system. 
Each RP is approaching the beam in small steps until a large signal increase is seen in the Beam 
Loss monitors directly downstream indicating that the RP is out of the shadow of the collimators 
and thus at the same number of (Jheam from the beam center as the collimators. The beam 
center position obtained from the collimator-based alignment is verified and corrected vertically 
and horizontally using the vertical position distribution from all particles in the horizontal pot and 
using the horizontal position distribution from a pure sample of elastic candidates, respectively. 
As a result of these collimator-based alignment exercises, the LHC operators has gained confidence 
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FIG. 3: The principle of the conimator-based RP ahgnment. The coUimators scrape the beam (left) creating 
a beam with sharp edges (second left). The top RP approaches the beam in small steps (middle) until large 
signal seen in beam loss monitor (second right). The procedure is repeated for the bottom RP (right). 

on the absolute positions of the RPs to allow the vertical (horizontal) RPs to approach closer to 
the beam in steps, taking data at 30 (30), 25 (30), 20 (25) and 18 (20) times the beam size, ahearm 
in standard LHC runs at y/s = 3.5 TeV and (3* = 3.5 m with nominal bunch intensity. Most (~ 
3.9 pb~^) of the integrated luminosity was taken with the vertical and horizontal RPs at IScr^eam 
and 20cj;,eam! respectively. In addition, the RPs have taken data with vertical and horizontal RPs 
at 1(Tbeam and IGcjfceam) respectively, during a collimator-based alignment exercise and during a 
dedicated few hours TOTEM run at the end of the data taking of 2010 corresponding to a total 
integrated luminosity of 9.5 nb^^. In the latter, a combined RP and T2 data taking was performed 
on a specially introduced bunch with a reduced number of protons {'^ 10^") for creating collisions 
without pileup for precise studies of forward charged particle multiplicity and diffractive processes. 

III. PHYSICS OF THE 2010 RUNS AND POTENTIAL OF THE 2011 RUNS 

The diffractive studies in TOTEM covers several diffractive processes: elastic scattering pp 
pp, single pp ^ p + X, double pp ^ Y + X and central diffraction pp ^ p + X + p, where + 
indicates a rapidity gap, an r/-range without final state particles, and X and Y hadronic systems. 

The 2010 LHC optics conditions of ^/s = 3.5 TeV and (3* = 3.5 m allow both protons in an elastic 
pp collision be detected by sector 45 top(bottom) and sector 56 bottom(top) RP combinations if 
\ty\ > ~0.4 and \ty\ > ~2.5 GeV^ in runs with vertical RPs at 7aheam and ISaheam, respectively. 
Diffractive protons can be seen in the 2010 LHC optics conditions either in the horizontal RPs (and 
possibly also in either top or bottom RPs if \ty\ is sufficiently large) if their ^ > ~2 % or only the 
vertical RPs with similar \ty\ requirements as elastic protons (if ^ < ~2 %). By selecting protons 
reconstructed in the horizontal RPs or only in the vertical RPs, one can enforce certain kinematics 
on the diffractive final state. As a demonstration, high and low mass single diffractive (SD) events 
can be selected by requiring the proton to be reconstructed in the horizontal RP or only in the 
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FIG. 4: Top: A high mass single diffraction candidate from data taken in coUisions of bunches with reduced 
number of protons October 2010. Bottom: A low mass single diffraction candidate from the same run. 

vertical RP, respectively, as shown in Fig. [H This kinematical selection can then be confirmed by 
the charged particles reconstructed in the T2 measuring charged particles 5.3 < \ri\ < 6.5 on both 
sides of IP5 with high mass SD having charged particles in both sides and low mass SD having 
charged particles only on the opposite side with respect to the proton (see Fig. H]). The expected 
physics results of the full analysis of the data taken by TOTEM in 2010 at ^/s = 7 TeV are: 

• a determination of the elastic scattering cross-section for ~0.4 < \t\ < ~5 GeV^; 

• a study of the kinematics of central diffractive events in the kinematically accessible range 
i.e. both p's with large .^'s, one p with large |i| and other with large or both p's with large 
|t|'s and a determination of the mass distribution of high mass central diffractive events; 

• a study of the ^ and rapidity gap correlation for single diffractive events and a determination 
of the mass distribution for high mass single diffractive events combining RP and T2 data; 
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• a determination of the forward charged particle multiphcity in the 5.3 < \r]\ < 6.5 region; 

• a study of the charged particle multiplicity correlation over a large rapidity range i.e. the 
forward-backward charged particle multiplicity correlations in inelastic pp collisions. 

The analysis of the data and the estimation of the systematic errors affecting these measurements 
is still in progress. The RP measurements on one hand are especially affected by uncertainties 
related to the LHC optics, the RP alignment and the trigger efficiencies and the T2 measurements 
on the other hand by uncertainties related to production of secondaries from interactions with 
material in front of the T2. Of the physics topics listed above, only the determination of the elastic 
cross-section in the high end of the accessible |t|-spectrum will be statistics limited. Therefore, in 
early 2011, the collimator-based RP alignment should be repeated with the new optics in nominal 
beam conditions to enable constant data taking with the RPs at 18cjf,eam or closer in normal LHC 
runs to improve significantly the statistics of elastic scattering at the highest accessible |t|-values. 

The installation of the Tl and the RP147 in the winter shutdown of 2010/11 and the efforts 
of the LHC machine to prepare for high /3* running in 2011 enables TOTEM to start pursuing 
its full physics programme during 2011. Having the Tl increases the charged particle detection 
range considerably from 5.3 < \rf\ < 6.5 to 3.1 < \r)\ < 6.5, enabling the inelastic events to be 
classified into non-diffractive minimum bias, single and double diffractive events with high purity 
and efficiency. A low P* run with collisions of bunches with reduced number of protons will allow 
the determination of the inelastic rate to a few % and individual process rates at 5-10 %, and to 
extend the forward charged particle studies to the range 3.1 < \r]\ < 6.5. At /3* = 90 m, all protons 
with \t\ > ~0.1 GeV^ will be detected by the RP220s irrespective of their ^ if the vertical RPs are 
approached to labeam- Short runs at /3* = 90 m opens up the possibility for following physics: 

• a determination of the total pp cross-section atot and luminosity at 5 % and 7 % level, 
respectively, using the luminosity-independent method. In addition, in conjunction with the 
atot measurement, the elastic and inelastic cross sections will be determined at 5 % level; 

• a comprehensive study of soft singel and central diffraction at any diffractive mass value. 

IV. SUMMARY AND CONCLUSIONS 

The TOTEM experiment has successfully commissioned its Roman Pot (RP) detectors at 220 
m and its forward inelastic telescope T2. In 2010, TOTEM has also taken data with the vertical 
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(horizontal) RP detectors at a distance of 7ai,eam i^^'^beam) in special and at a distance of IScTfteam 
{20abeam) in nominal LHC runs allowing TOTEM to complete in its first physics measurement, the 
measurement of the elastic pp cross-section at s/s = 7 TeV for the 0.4 < |t| < 5 GeV^ range, in the 
upcoming months. In addition, TOTEM will in the near future produce results on high mass single 
and central diffraction as well as the forward charged particle multiplicity and forward- backward 
charged particle multiplicity correlation for in the 5.3 < jr/j < 6.5 range. 

With the installation of the Tl telescope in the 2010/2011 winter shutdown and the preparation 
of the /3* = 90 m optics from the LHC machine side, TOTEM will be allowed to start pursuing its 
full physics programme in 2011 by taking data in a combination of nominal LHC runs, dedicated 
high P* TOTEM runs and special low P* runs with bunches with reduced number of protons. /3* = 
90 m runs will allow a first measurement of the total pp cross-section at 5 % level and in conjunction 
measurements of the luminosity, elastic and inelastic cross-sections at 5-7 % level. In addition, P* 
= 90 m runs will allow a comprehensive study of single and central diffractivc processes at any mass 
as well as a determination of the cross-section of non-diffractivc minimum bias, single and double 
diffraction at the 5-10 % level. The forward charged particle multiplicity in the 3.1< \r]\ <6.5 
region and corresponding forward-backward multiplicity correlation will be studied either in a P* 
= 90 m run or in a special low P* run with bunches with reduced number of protons. 
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LHCf and the connection with Cosmic Ray Physics 
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The LHCf experiment has been designed to precisely measure the 7 and n energy spectra in 
the very forward region at LHC; these measurements are important to calibrate the Monte 
Carlo models widely used in the High Energy Cosmic Ray induced air showers analysis, 
allowing a better understanding of the experimental systematics. LHCf has started the data 
taking at the end of 2009, at 900 GeV center of mass energy, and has later on continued the 
data taking during the hrst part of 2010, when LHC center of mass energy was increased 
up to 7 TeV. Preliminary results of 900 GeV and 7 TeV data taking periods are reported in 
this paper. 







I. LHCF PHYSICS CASE 

The LHCf experiment is devoted to the precise measurements of the energy and transverse mo- 
mentum spectra of 7 and neutrons produced in the very forward region proton-proton interactions 
at the LHC collider, in the |r/| > 8.3 pseudo-rapidity region. The detector consists of a double 
arm (Arml and Arm2) - double tower sampling and imaging calorimeter, placed at it 140 m from 
ATLAS interaction point (IPl) inside the zero-degree neutral absorbers (Target Neutral Absorber, 
TAN). Charged particles from the IP are swept away by the inner beam separation dipole before 
reaching the TAN, so that only photons mainly from vr^ decays, neutrons and neutral kaons reach 
the LHCf calorimeters. 

Each calorimeter tower is made of 16 layers of plastic scintillators interleaved by tungsten layers 
as converter, complemented by a set of four X-Y position sensitive layers which provide incident 
shower positions, in order to obtain the transverse momentum of the incident primary and to cor- 
rect for the effect of leakage from the edges of the calorimeters. 

While the two calorimeters are identical for the calorimetric structure, they slightly differ for the 
geometrical arrangement of the two towers and for the position sensitive layers made by 1 mm^ scin- 
tillating fibers in one calorimeter and silicon micro-strip layers in the other. A detailed description 
of the LHCf detector and it's performances can be found in Qj. 



'Electronic address: Iadriani@fi.infti.it1 
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The LHCf experiment differs from the other LHC experiment not only in dimensions but also 
for the main physics motivation, which for LHCf is strictly connected to astroparticle physics. 
The goal of the experiment is indeed to measure neutral particle spectra to calibrate Monte Carlo 
codes used in High Energy Cosmic Ray (HECR) physics. A good knowledge of nuclear interaction 
model of primary cosmic rays with earth's atmosphere is mandatory to better understand many 
properties of primary cosmic rays, like the energy spectrum and the composition, whose knowledge 
is finally strictly related to our capability to understand the origin of high energy phenomena in the 
Universe. Dedicated extensive air shower experiments are taking data since many years and have 
strongly contributed to our understanding of High and Ultra High Energy Cosmic Ray (UHECR) 
Physics. However, the results of these experiments are in some cases not fully in agreement and, 
in addition, the interpretation of their data in terms of primary cosmic ray properties is strongly 
affected by the knowledge of the nuclear interactions in the earth's atmosphere. This is true, for 
instance, for the interpretation of the behaviour of the energy spectrum in the UHE region, in 
particular the existence of events above the so called GZK cut-off, and the chemical composition 
of cosmic rays. Indeed, evidence of UHECR, above the GZK cut-off, has been reported for the 
first time by the AGASA experiment On the contrary, the results of the HiRes J] experiment 
and, more recently, the ones of the Pierre Auger Collaboration [s] are consistent with the existence 
of the cut-off. The disagreement among data would be reduced by adjusting the energy scales 
of the different experiments to account for systematic effects in the determination of the particle 
energy, that might be due to different detecting techniques. Similar considerations hold for the 
interpretation of cosmic ray composition since it is directly related to their primary sources. Ac- 
celerator experiments validating the interaction model chosen are hence essential. As a matter of 
fact air shower development is dominated by the forward products of the interaction between the 
primary particle and the atmosphere. The only available data on the production cross-section of 
neutral pions emitted in the very forward region have been obtained more than twenty years ago 
by the UA7 Collaboration [g] at the CERN SppS up to an energy of 10^^ eV and in a very narrow 
pseudo-rapidity range. The LHCf experiment at LHC has the unique opportunity to take data at 
energies ranging from ^/s = 0.9 TeV up to 14 TeV, thus extending significantly the energy range 
up to a region of great interest for high energy cosmic rays, the region between the "knee" and the 
GZK cut-off. 

Additionally, the possibility to install the detector in the region where the single LHC beam 
pipe splits in two allows us to cover the pseudo-rapidity range \ri\ > 8.3, where most of the energy 
flux produced in the p-p collisions is concentrated. Fig. [1] shows the distribution of the number 
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of particles (left side) and of the energy flux (right side) as function of the pseudo-rapidity; most 
of the final state particles in the p-p reaction are produced in the central region, but most of the 
energy flux is concentrated in the very forward region, 6 < \ri\ < 12, clearly demonstrating the 
advantage of a very forward detector in the energy flux measurement. 




FIG. 1: Distributions of the number of particles (left side) and of the energy flux (right side) as function 
of pseudo-rapidity in the p-p collisions at 14 TeV center of mass energy. Typical coverage of the various 
central and forward detectors are also shown. 



II. THE EXPERIMENTAL ASPECTS 

In order to calibrate the Monte Carlo codes used in HECR physics, the LHCf experiment should 
be able to have a detailed knowledge of the absolute energy scale. For this reason LHCf relies on 
a very precise reconstruction of the vr" mass, by reconstructing in the two towers the showers from 
the 2 7 from vr'^ decays. 

The performances of the detector have been careful measured using beam test data and well satisfy 
the design requirements l|. The measured position resolution in locating the shower center for 
particles above 100 GeV (which is the region of interest for LHCf) is about 200 fim for scintillating 
fibre layers (ARMl) and about 50 fim for the silicon micro-strip layers (ARM2). 
The energy resolution for the calorimeters is better than 4% at 200 GeV, as can be seen from 
Fig. [21 that shows the energy resolution for electromagnetic particles measured on beam test and 
expected from the Monte Carlo simulation, for two different photomultipliers high voltage setting 
(Low Gain and High Gain mode). This figure demonstrate the excellent performances of LHCf in 
the high energy electromagnetic particles reconstruction, despite the small transverse size of the 
towers (20 x 20 mm^ and 40 x 40 mm^). 
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FIG. 2: Energy resolution for electromagnetic particles measured on beam test (red markers) and expected 
from the Monte Carlo simulation (empty markers). The results are shown for two different setting of the 
photomultipliers high voltage, corresponding to Low Gain and High Gain mode. 

Fig. [3] shows LHCf expectations for the 7 and neutron energy spectra for few minutes exposure 
at 10^^ cm^^s^^ with 3.5+3.5 TeV center of mass energy p-p collisions. Depending on the nuclear 
interaction model used, the energy spectra change more or less significantly. As can be seen from 
this plot, the LHCf experiment will be able to disentangle different interaction models already at 
lower energy and with very low statistics, thus ensuring a calibration of cosmic ray Monte Carlo 
in an energy range wider than the one expected at the beginning of the project. 

Gamma \/s=7TeV Hadron\fs=7TeV 




500 1000 1500 2000 1000 2000 3000 4000 5000 

Energy(G Energy(GeV) 

FIG. 3: Expected energy spectrum for photons (left) and neutral hadrons (right) according to different 
interaction models at 3.5+3.5 TeV center of mass energy p-p collisions. 
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III. DATA TAKING AT 900 GEV 



At the end of November 2009 LHC has started to provide colhsions to the experiments at 900 
GeV center of mass energy. The LHCf experiment has taken data from December 6 till December 
15, accumulating about 6500 shower triggers in total on both arm of the calorimeter. The hit maps 
obtained in a typical run for Arml and Arm2 detectors are shown in Fig clearly showing the 
effect on the detector acceptance of the elliptical beam pipe close to the bending dipole magnet 
Dl in the region between interaction point and the LHCf location. A typical 7 event registered on 
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FIG. 4: Hit maps obtained in a typical LHCf run for tlie Arml (left side) and Arm2 (right side) detectors. 
The acceptance reduction due to the elliptical beam pipe located close to the bending dipole magnet Dl in 
the region between interaction point and the LHCf location is clearly visible. 

ARM2 detector is shown in Fig. [5l 

A preliminary analysis has been carried on to reconstruct 7 and hadron spectra. The particle 
identification has been achieved through the use of transition curve information. Results obtained 
in the two towers as well as in the two arms are consistent each other both for photons and hadrons 
as can be seen from Fig. [6l 

900 GeV data have been accumulated also in a second run during spring 2010 and whole analysis 
of the events collected in 2009 and 2010 is ongoing and almost ready for publication. 
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FIG. 5: A typical 7 event registered on ARM2 detector in 900 GeV collision data. The two upper panels 
show the longitudinal energy profile deposited on each tower of the calorimeter, while the two lower panels 
show the transverse energy X and Y profile deposited on each of the four silicon layers. 

IV. DATA TAKING AT 7 TEV 

At the beginning of 2010 the LHC beam energy was increased up to 3.5 TeV, allowing LHCf 
to take data at 7 TeV center of mass energy. The geometrical configuration of the detectors and 
the kinematic of the decays allow LHCf to reconstruct vr" events, by measuring the two 7 of the 
decay in the two separate towers. In this way the absolute energy scale calibration can be cross 
checked by looking at the kinematically reconstructed vr*^ invariant mass. A typical vr*^ candidate 
event collected in ARM2 detector is shown in Fig. [71 

Data analysis for the 7 TeV running is still ongoing. As preliminary result. Fig. [8] shows the 
measured vr*^ invariant mass and energy spectra obtained for Arml and Arm2. The invariant mass 
distributions demonstrate the excellent performances in the vr" reconstruction (Am/m ~ 5% for 
Arml and ~ 2.3% for Arm2), even for the extremely high energy tt*^ (up to 3.5 TeV). 

V. FUTURE ACTIVITIES 

LHCf has been originally designed to take data at high energy (14 TeV) and low luminosity (L~ 
10^^ cm^^s^^), in the early LHC operation phase. However, the 2008 LHC incident significantly 
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FIG. 6: Comparison of reconstructed 7 and hadron spectra in small and large tower of ARM2 detector 
(upper plot) and in ARMl and ARM2 detector (bottom plot) for 900 GeV p-p data. Spectra are normalized 
to take into account different geometrical acceptance. 

changed the running plans, with an extended initial running period at reduced energy (7 TeV) and 
quite high luminosity (> 10^^ cm~^s~^). The plastic scintillators used for the energy measurement 
are intrinsically radiation weak, and can not sustain for long periods the luminosity foreseen during 
the 2010 running. Fig. [9] shows the results of the measurements that have been done on the plastic 
scintillators with the heavy ion beam from a synchrotron at HIMAC of NIRS (National Institute of 
Radiological Science, Japan) and with 7 rays at the ^"^Co Radiation Facility of Nagoya University. 
The light output reduction for the EJ260 plastic scintillators was measured as function of the 
integrated dose, showing a significant reduction in the light yield for doses greater than few tens of 
Grays. In the normal LHC running conditions at 7 TeV, a 10 Gy dose is gathered for an integrated 
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FIG. 7: A TT*^ candidate event registered on ARM2 detector in 7 TeV collision data. The plots follow the 
same convention as in Fig. [5] 
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FIG. 8: Measured invariant mass and energy spectra for 7r° at 7 TeV center of mass energy, for Arml (top 
plots) and Arm2 (bottom plots). The invariant mass resolution is of the order of 5% and 2.3% respectively. 
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luminosity of the order of 20 nb ^, corresponding to few days of data taking at > 10^^ cm ^s ^ 
luminosity. 
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FIG. 9: The variation of the scintillator light outputs for various types of plastic scintillators, including 
EJ260, as function of the dose integrated by using heavy ions and 7 rays. 

The detectors were hence removed at the end of July 2010, during a 3 days LHC technical stop, 
after a ~ 10 Gy accumulated dose. They will be tested at the H4 SPS beam test facility during 
Autumn 2010, to confirm the absolute energy scale. Later on, during 2011, the detectors will be 
upgraded, by replacing plastic scintillators with radiation hard GSO crystals, able to sustain doses 
of ~ 10^ — 10'' Gy without significant damages. In this way, LHCf will be able to be installed again 
in the tunnel in 2013, when the LHC energy will be increased up to 14 TeV. 



VI. CONCLUSIONS 



The LHCf experiment has taken data at LHC both at 900 GeV as well at 7 TeV center of mass 
energy, before being removed at the end of July 2010. It will be upgraded with radiation hard 
scintillators, and re-installed in the LHC tunnel in 2013, when the machine will provide 7+7 TeV 
p-p collisions. 
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Thanks to the excellent detector performance for the reconstruction of photon, neutral meson 
and neutron spectra at different energies from 900 GeV up to 14 TeV p-p runs, LHCf will allow 
to calibrate air shower Monte Carlo codes covering the most interesting energy range for HECR 
Physics, thus providing invaluable input to our understanding of high energy phenomena in the 
Universe. 
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Mueller Navelet jets, jet gap jets and anomalous WWj'j couplings in 7-induced 

processes at the LHC 

C. Royon^' 

^ IRFU/Service de physique des particules, 
CEA/Saclay, 91191 Gif-sur-Yvette cedex, France 

We describe two different important measurements to be performed at the LHC. The 
Mueller Navelet jet and jet gap jet cross section represent a test of BFKL dynamics and we 
perform a NLL calculation of these processes and compare it with recent Tevatron measure- 
ments. The study of the WW'yy couplings at the LHC using the forward detectors proposed 
in the ATLAS Forward Physics project as an example allows to probe higgsless and extradi- 
mension models via anomalous quartic couplings since the reach is improved by four orders 
of magnitude with respect to the LEP results.. 



I. MUELLER NAVELET JETS AT THE LHC 

In this section, we give the BFKL NLL cross section calculation for Mueller Navelet processes 
at the Tevatron and the LHC. Since the starting point of this study was the description of forward 
jet production at HERA, we start by describing briefly these processes. 

A. Forward jets at HERA 



Following the successful BFK 



Leading Order (LO) at HERA 



m 



parametrisation of the forward-jet cross-section da/dx at 



it is possible to perform a similar study using Next-to-leading 
(NLL) resummed BFKL kernels. Forward jets at HERA are an ideal observable to look for BFKL 
resummation effects. The interval in rapidity between the scattered lepton and the jet in the 
forward region is large, and when the photon virtuality is close to the transverse jet momentum 
kx, the DDLAP cross section is small because of the fcy ordering of the emitted gluons. In this 
short report, we will only discuss the phenomelogical aspects and all detailed calculations can be 

1 1 rf 

found in Ref. [J] for forward jets at HERA and in Ref. [5] for Mueller Navelet jets at the Tevatron 
and the LHC. 
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FIG. 1: Comparison between the HI da/dx measurement with predictions for BFKL-LL, BFKL-NLL (S3 
and S4 schemes) and DGLAP NLO calculations (see text). S4, S3 and LL BFKL cannot be distinguished 
on that figure. 



B. BFKL NLL formalism 



The BFKL NLL [6| longitudinal transverse cross section reads: 



d(j> 



dxjdkj, 

where xj is the proton momentum fraction carried by the forward jet, Xeff is the effective BFKL 
NLL kernel and the (ps are the transverse and longitunal impact factors taken at LL. The effective 
kernel Xe//(7)«) is defined from the NLL kernel Xnll{i-,^) by solving the implicit equation 
numerically 



Xeffin, a) = XNLL b, OL Xeffil, «)] 



(2) 



The integration over 7 in Eq. [T] is performed numerically. It is possible to fit directly da/dx 
measured by the HI collaboration using this formalism with one single parameter, the normali- 
sation. The values of xnll ai"e taken at NLL 6^ using different resummation schemes to remove 
spurious singularities defined as S3 and S4 [7]. Contrary to LL BFKL, it is worth noticing that the 
coupling constant as is taken using the renormalisation group equations, the only free parameter 
in the fit being the normalisation. 

To compute da/dx in the experimental bins, we need to integrate the differential cross section on 
the bin size in Q^, xj (the momentum fraction of the proton carried by the forward jet), /ct , while 
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d a/dx dp^2 d - HI DATA 




FIG. 2: Comparison between the HI measurement of the triple differential cross section with predictions for 
BFKL-LL, BFKL-NLL and DGLAP NLO calculations (see text). 

taking into account the experimental cuts. To simplify the numerical calculation, we perform the 
integration on the bin using the variables where the cross section does not change rapidly, namely 
k^/Q'^, logl/xj, and l/Q^. Experimental cuts are treated directly at the integral level (the cut 
on 0.5 < kj-ZQ"^ < 5 for instance) or using a toy Monte Carlo. More detail can be found about the 
fitting procedure in Appendix A of Ref. [3(]. 

The NLL fits fl] can nicely describe the HI data |3] for the S4 and S3 schemes ^-4\ (x^ = 0.48/5 
and = 1.15/5 respectively per degree of freedom with statistical and systematic errors added in 
quadrature) . The curve using a LL fit is indistinguishable in Fig. [1] from the result of the BFKL- 
NLL fit. The DGLAP NLO calculation fails to describe the HI data at lowest x (see Fig. [T]). We 
also checked the effect of changing the scale in the exponential of Eq. [1] from kxQ to Ik^Q or 
kxQ/^ which leads to a difference of 20% on the cross section while changing the scale to k^ or 
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modifies the result by less than 5% which is due to the cut on 0.5 < K'^jQ'^ < 5. Implementing 
the higher-order corrections in the impact factor due to exact gluon dynamics in the 7* — t- qq 
transition changes the result by less than 3%. ^ 

The HI collaboration also measured the forward jet triple differential cross section [8] and the 
results are given in Fig. [21 We keep the same normalisation coming from the fit to da / dx to predict 
the triple differential cross section. The BFKL LL formalism leads to a good description of the 
data when r = k^/Q'^ is close to 1 and deviates from the data when r is further away from 1. This 
effect is expected since DGLAP radiation effects are supposed to occur when the ratio between the 
jet kx and the virtual photon are further away from 1. The BFKL NLL calculation including 
the evolution via the renormalisation group equation leads to a good description of the HI 
data on the full range. We note that the higher order corrections are small when r ~ 1, when the 
BFKL effects are supposed to dominate. By contrast, they are significant as expected when r is 
different from one, ie when DGLAP evolution becomes relevant. We notice that the DGLAP NLO 
calculation fails to describe the data when r ~ 1, or in the region where BFKL resummation effects 
are expected to appear. 

In addition, we checked the dependence of our results on the scale taken in the exponential of 
Eq. [TJ The effect is a change of the cross section of about 20% at low pT increasing to 70% at 
highest pt- Taking the correct gluon kinematics in the impact factor lead as expected to a better 
description of the data at high px . 



C. Mueller Navelet jets at the Tevatron and the LHC 



Mueller Navelet jets are ideal processes to study BFKL resummation effects [10(]. Two jets with 
a large interval in rapidity and with similar tranverse momenta are considered. A typical observable 
to look for BFKL effects is the measurement of the azimuthal correlations between both jets. The 
DGLAP prediction is that this distribution should peak towards vr - ie jets are back-to-back- 
whereas multi-gluon emission via the BFKL mechanism leads to a smoother distribution. The 
relevant variables to look for azimuthal correlations are the following: 

^r] = yi- y2 

y = {yi + 2/2)72 



Q = \/kik2 
R = k2/ki 
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FIG. 3: The Mueller-Navelet jet A"I> distribution for LHC kinematics in the BFKL framework at LL (upper 
plots) and NLL-S4 (lower plots) accuracy for A77 = 6, 8, 10. 
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FIG. 4: Azimuthal correlations between jets with A77 ^=6, 8, 10 and 11 and pr > 5 GeV in the CDF 
acceptance. This measurement will represent a clear test of the BFKL regime. 



where yi^2 and ki^2 are respectively the jet rapidities and transverse momenta. The azimuthal 
correlation for BFKL reads: 



27r 
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where in the NLL BFKL framework, 



^as{QVR)as{Q^R) {^j"^ dyxif,ff{xi,Q^/R)x2feff{x2,Q^R) 



1/2+00 ^ 

1/2-00 2i7r 

and Xeff is the effective resummed kernel. Computing the different CTp at NLL for the resummation 
schemes S3 and S4 aUowed us to compute the azimuthal correlations at NLL. As expected, the A$ 
dependence is less flat than for BFKL LL and is closer to the DGLAP behaviour j^. In Fig. [3l 
we display the observable 1/ adcr / d/\^ as a function of A<I>, for LHC kinematics. The results are 
displayed for different values of A?? and at both LL and NLL accuracy using the S4 resummation 
scheme. In general, the A$ spectra are peaked around A$ = 0, which is indicative of jet emissions 
occuring back-to-back. In addition the A<I> distribution flattens with increasing Ar] = yi—y2. Note 
the change of scale on the vertical axis which indicates the magnitude of the NLL corrections with 
respect to the LL-BFKL results. The NLL corrections slow down the azimuthal angle decorrelations 
for both increasing An and R deviating from 1. We also studied the R dependence of our prediction 
which is quite weak [s] and the scale dependence of our results by modifying the scale to either 
or 2Q^ and the effect on the azimuthal distribution is of the order of 20%. The effect of the 

n 

energy conservation in the BFKL equation [5|] is large when R goes away from 1. The effect is to 
reduce the effective value of Ai] between the jets and thus the decorrelation effect. However, it is 
worth noticing that this effect is negligible when R is close to 1 where this measurement will be 
performed. 

A measurement of the cross-section da^^'^^'^ /dArjdRdA^ at the Tevatron (Run 2) or the LHC 
will allow for a detailed study of the BFKL QCD dynamics since the DGLAP evolution leads to 
much less jet angular decorrelation (jets are back-to-back when R is close to 1). In particular, 
measurements with values of Arj reaching 8 or 10 will be of great interest, as these could allow to 
distinguish between BFKL and DGLAP resummation effects and would provide important tests 
for the relevance of the BFKL formalism. 

To illustrate this result, we give in Fig. H] the azimuthal correlation in the CDF acceptance. 
The CDF collaboration installed the mini-Plugs calorimeters aiming for rapidity gap selections in 
the very forward regions and these detectors can be used to tag very forward jets. A measurement 
of jet pt with these detectors would not be possible but their azimuthal segmentation allows a (p 
measurement. In Fig. SI we display the jet azimuthal correlations for jets with a pT > 5 GeV and 
Ar] =6, 8, 10 and 11. For Ar] =11, we notice that the distribution is quite flat, which would be a 
clear test of the BFKL prediction. 
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II. JET GAP JETS AT THE TEVATRON AND THE LHC 

In this section, we describe another possible measurement which can probe BFKL resummation 
effects and we compare our predictions with existing DO and CDF measurements [ll]. 

A. BFKL NLL formalism 

The production cross section of two jets with a gap in rapidity between them reads 



where ^/s is the total energy of the collision, Et the transverse momentum of the two jets, xi and 
X2 their longitudinal fraction of momentum with respect to the incident hadrons, S the survival 

n 

probability, and / the effective parton density functions [11]. The rapidity gap between the two 
jets is A7/ = ln(xiX2s/py). 

The cross section is given by 

^^^ = ^|A,A„,£|)r (4) 

in terms of the gg — )• gg scattering amplitude ^(Ar/,p|,). 

In the following, we consider the high energy limit in which the rapidity gap A?? is assumed to 
be very large. The BFKL framework allows to compute the gg — t- gg amplitude in this regime, and 
the result is known up to NLL accuracy 

A(/\ p2^ leiVcTTol ^ r d7 [p'-(7-l/2)']exp{a(i?|)xe//[2p,7,«(^j-)]A^} , . 
^ CpEl ^£^^J 2^n [(7 - 1/2)^ - - 1/2)2] [(7 - 1/2)2 - (p + 1/2)2] 

with the complex integral running along the imaginary axis from 1/2 — ioo to 1/2 + zoo, and with 
only even conformal spins contributing to the sum, and a = osNc/tt the running coupling. 

Let us give some more details on formula [5j The NLL-BFKL effects are phenomenologically 
taken into account by the effective kernels Xeff{p,1,o). The NLL kernels obey a consistency 
condition which allows to reformulate the problem in terms of Xef fil^ot)- The effective kernel 
Xe//(7;Q;) is obtained from the NLL kernel Xnll{1i^) by solving the implicit equation Xeff = 
XMLiin^oi Xeff) as a solution of the consistency condition as it was also performed for forward 
jets. 

In this study, we performed a parametrised distribution of da^^'^^^ /dE^ so that it can be easily 
implemented in the Herwig Monte Carlo {l^ since performing the integral over 7 in particular would 
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FIG. 5: Comparisons between the DO measurements of the jet-gap-jet event ratio with the NLL- and LL- 
BFKL calculations. The NLL calculation is in fair agreement with the data. The LL calculation leads to a 
worse description of the data. 



be too much time consuming in a Monte Carlo. The implementation of the BFKL cross section 
in a Monte Carlo is absolutely necessary to make a direct comparison with data. Namely, the 
measurements are sensititive to the jet size (for instance, experimentally the gap size is different 
from the rapidity interval between the jets which is not the case by definition in the analytic 
calculation) . 



B. Comparison with DO and CDF measurements 

Let us first notice that the sum over all conformal spins is absolutely necessary. Considering 
only p = in the sum of Equation [5] leads to a wrong normalisation and a wrong jet Et dependence, 
and the effect is more pronounced as Ar/ diminishes. 

The DO collaboration measured the jet gap jet cross section ratio with respect to the total 

dijet cross section, requesting for a gap between -1 and 1 in rapidity, as a function of the second 

leading jet Et, and Arj between the two leading jets for two different low and high Et samples 

(15< Et <20 GeV and Et >30 GeV). To compare with theory, we compute the following quantity 

_ BFKL NLL HERWIG LO QCD 
~ Dijet Herwig NLO QCD 
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FIG. 6: Comparisons between the CDF measurements of the jet-gap-jet event ratio with the NLL- and 
LL-BFKL calculations. The NLL calculation is in fair agreement with the data. The LL calculation leads 
to a worse description of the data. 



in order to take into account the NLO corrections on the dijet cross sections, where BFKL NLL 
HERWIG and Dijet Herwig denote the BFKL NLL and the dijet cross section implemented in 
HERWIG. The NLO QCD cross section was computed using the NLO Jet -|— I- program [l^. 

The comparison with DO data is shown in Fig. 5. We find a good agreement between 
the data and the BFKL calculation. It is worth noticing that the BFKL NLL calculation leads 
to a better result than the BFKL LL one (note that the best description of data is given by the 
BFKL LL formalism for p = but it does not make sense theoretically to neglect the higher spin 
components and this comparison is only made to compare with previous LL BFKL calculations). 

The comparison with the CDF data [l^ as a function of the average jet Ex and the difference 
in rapidity between the two jets is shown in Fig. 6, and the conclusion remains the same: the 
BFKL NLL formalism leads to a better description than the BFKL LL one. 



Predictions for the LHC 



Using the same formalism, and assuming a survival probability of 0.03 at the LHC, it is possible 
to predict the jet gap jet cross section at the LHC. While both LL and NLL BFKL formalisms 
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FIG. 7: Ratio of the jet gap jet to the inchisive jet cross sections at the LHC as a function of jet pr and At]. 



lead to a weak jet Et or Ar] dependence, the normalisation is found to be quite different (see Fig. 
7) leading to higher cross section for the BFKL NLL formalism. 
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III. QUARTIC ANOMALOUS COUPLINGS AT THE LHC 

In the third part of this report, we discuss a completely different topic, namely the possibility 
to probe anomalous quartic couplings between photons and W or Z bosons at the LHC with an 
unprecedent precision using forward detectors to be installed in CMS and ATLAS experiments [3]. 
In the Standard Model (SM) of particle physics, the couplings of fermions and gauge bosons are 
constrained by the gauge symmetries of the Lagrangian. The measurement of W and Z boson pair 
productions via the exchange of two photons allows to provide directly stringent tests of one of 
the most important and least understood mechanism in particle physics, namely the electroweak 
symmetry breaking [l^. The non-abelian gauge nature of the SM predicts the existence of quartic 
couplings WWjj between the W bosons and the photons which can be probed directly at the 
Large Hadron Collider (LHC) at CERN. The quartic coupling to the Z boson ZZ^^ is not present 
in the SM. Quartic anomalous couplings between the photon and the Z or W bosons are specially 
expected to occur in higgsless or extradimension models jl?]. 

A. Photon exchange processes in the SM 
The process that we intend to study is the W pair production shown in Fig. 8 induced by the 



exchange of two photons 



15 



18l |. It is a pure QED process in which the decay products of the W 
bosons are measured in the central detector and the scattered protons leave intact in the beam pipe 
at very small angles, contrary to inelastic collisions. Since there is no proton remnant the process 
is purely exclusive; only W decay products populate the central detector, and the intact protons 
can be detected in dedicated detectors located along the beam line far away from the interaction 
point. 

The cross section of the pp — t- pWWp process which proceeds through two-photon exchange is 
calculated as a convolution of the two-photon luminosity and the total cross section 77 — )• WW. 
The total two-photon cross section is 95.6 fb. 

All considered processes (signal and background) were produced using the Forward Physics 
Monte Carlo [20,] (FPMC) generator. The aim of FPMC is to produce different kinds of processes 
such as inclusive and exclusive diffractio n, p hoton-exchange processes. FPMC was interfaced to as 



fast simulation of the ATLAS detector 2l|]. To reduce the amount of considered background, we 
only use leptonic (electrons and muons) decays of Z and W bosons. The following backgrounds 
were considered: 77 — )■ — two-photon dilepton production, DPE— )• — - dilepton production 
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FIG. 8: Sketch diagram showing the two-photon production of a central system. 
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FIG. 9: Distribution of the transverse momentum of the leading lepton for signal and background after the 
cut on W, ^T, and Acj) between the two leptons. 

through double pomeron exchange, DPE— )• W/'+VF^ — >• lli'i' — diboson production through double 
pomeron exchange. 

After simple cuts to select exclusive W pairs decaying into leptons, such as a cut on the proton 
momentum loss of the proton (0.0015 < ^ < 0.15) — we assume the protons to be tagged in 
the ATLAS Forward Physics detectors ^] — , on the transverse momentum of the leading and 
second leading leptons at 25 and 10 GeV respectively, on > 20 GeV, Acj) > 2.7 between leading 
leptons, and 160 <W< 500 GeV, the diffractive mass reconstructed using the forward detectors, 
the background is found to be less than 1.7 event for 30 fb^^ for a SM signal of 51 events. In this 
channel, a 5 a discovery of the Standard Model pp — )• pWWp process is possible after 5 fb^^. 
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B. Quartic anomalous couplings 



The parameterization of the quartic couphngs based on [22] is adopted. We concentrate on the 
lowest order dimension operators which have the correct Lorentz invariant structure and obey the 
SU{2)c custodial symmetry in order to fulfill the stringent experimental bound on the p parameter. 
The lowest order interaction Lagrangians which involve two photons are dim-6 operators. The 
following expression for the effective quartic Lagrangian is used 

_p2 W 2 Z 

where oq, ac are the parametrized new coupling constants and the new scale A is introduced so 
that the Lagrangian density has the correct dimension four and is interpreted as the typical mass 
scale of new physics. In the above formula, we allowed the W and Z parts of the Lagrangian to 
have specific couplings, i.e. — )• (flo^, ) ^^'^ similarly ac — )• {oq , 

The WW and ZZ two-photon cross sections rise quickly at high energies when any of the 
anomalous parameters are non-zero. The cross section rise has to be regulated by a form factor 
which vanishes in the high energy limit to construct a realistic physical model of the BSM theory. 
We therefore modify the couplings by form factors that have the desired behavior, i.e. they modify 
the coupling at small energies only slightly but suppress it when the center-of-mass energy W^^ 
increases. The form of the form factor that we consider is the following 

(8) 



(1 + W^^/A^)^ 
where n=2, and A ~2 TeV. 

The cuts to select quartic anomalous gauge coupling WW events are similar as the ones we 
mentioned in the previous section, namely 0.0015 < ^ < 0.15 for the tagged protons, > 20 GeV, 
A(f) < 3.13 between the two leptons. In addition, a cut on the pT of the leading lepton px > 160 
GeV and on the diffractive mass W > 800 GeV are requested since anomalous coupling events 
appear at high mass. Fig 9 displays the px distribution of the leading lepton for signal and the 
different considered backgrounds. After these requirements, we expect about 0.7 background events 
for an expected signal of 17 events if the anomalous coupling is about four order of magnitude lower 
than the present LEP limit (log^/A^I = 5.4 10~^) for a luminosity of 30 fb^^. The strategy to 
select anomalous coupling ZZ events is analogous and the presence of three leptons or two like sign 
leptons are requested. Table 1 gives the reach on anomalous couplings at the LHC for a luminosity 
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Couplings 


OPAL limits 
[GeV-2] 


Sensitivity 

5(7 


@ £ = 30 (200) fb-i 
95% CL 




[-0.020, 0.020] 


5.4 10-6 
(2.7 10-6) 


2.6 10-6 
(1.4 10-6) 




[-0.052, 0.037] 


2.0 10-5 
(9.6 10-6) 


9.4 10-6 
(5.2 10-6) 


af/A2 


[-0.007, 0.023] 


1.4 10-5 
(5.5 10-6) 


6.4 10-6 
(2.5 10-6) 


ag/A2 


[-0.029, 0.029] 


5.2 10-5 
(2.0 10-5) 


2.4 10-5 
(9.2 10-6) 



TABLE I: Reach on anomalous couplings obtained in 7 induced processes after tagging the protons in the 
final state in the ATLAS Forward Physics detectors compared to the present OPAL limits. The 5a discovery 
and 95% C.L. limits are given for a luminosity of 30 and 200 ft)-^ 



231 ] ■ We note that we can gain almost 



of 30 and 200 fb-^ compared to the present OPAL limits 
four orders of magnitude in the sensitivity to anomalous quartic gauge couplings compared to LEP 
experiments, and it is possible to reach the values expected in Higgsless or extra-dimension models 
which are of the order of 5 10"^. The tagging of the protons using the ATLAS Forward Physics 
detectors is the only method at present to test such small values of quartic anomalous couplings and 
thus to probe the higgsless models in a clean way. The reach on anomalous triple gauge couplings 
is much less improved at the LHC compared to LEP experiments {2^]. 

To conclude, the ATLAS Forward Physics program (and the CMS one) will allow to study 
Higgsless models with an unprecedent precision as well as to probe the Higgs boson by allowing 
its mass and spin measurements ^25i] using the forward detectors proposed for installation at 220 
and 420 m in ATLAS and CMS. 
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Experimental results on diffraction at CDF 

Michele Gallinaro* ^ 

^ The Rockefeller University (USA) and LIP Lisbon, Portugal 

Diffractive events are studied by means of identification of one or more rapidity gaps 
and/or a leading antiproton. Measurements of soft and hard diffractive processes have been 
performed at the Tevatron pp coUider and presented. We report on the diffractive structure 
function obtained from dijet production in the range < < 10, 000 GeV^, and on the \t\ 
distribution in the region < |f| < 1 GeV^ for both soft and hard diffractive events up to 
w 4,500 GeV^. Results on single diffractive W/Z production, forward jets, and central 
exclusive production of both dijets and Z-bosons are also presented. 



I. INTRODUCTION 

Diffraction can be described as an exchange of a combination of quarks and gluons carrying 
the quantum numbers of the vacuum jl]. As no radiation is expected from such an exchange, 
diffractive processes are characterized by the presence of large rapidity regions not filled with 
particles ("rapidity gaps"). 

At the Fermilab Tevatron collider, proton-antiproton collisions have been used to study diffrac- 
tive interactions in Run I (1992-1996) at an energy of ^/s = 1.8 TeV and continued in Run II 
(2003-present) with new and upgraded detectors at ^/s = 1.96 TeV. The goal of the CDF exper- 
imental program at the Tevatron is to provide results help decipher the QCD nature of hadronic 
diffractive interactions, and to measure exclusive production rates which could be used to es- 
tablish the benchmark for exclusive Higgs production at the Large Hadron Collider (LHC). The 
study of diffraction has been performed by tagging events either with a rapidity gap or with a 
leading hadron. The experimental apparatus includes a set of forward detectors [z] that extend 
the rapidity coverage to the forward region. The Miniplug (MP) calorimeters cover the region 
3.5 < \r]\ < 5.1; the Beam Shower Counters (BSC) surround the beam-pipe at various locations and 
detect particles in the region 5.4 < |r/| < 7.4; the Roman Pot spectrometer (RPS) tags the leading 
hadron scattered from the interaction point after losing a fractional momentum approximately in 
the range 0.03 < ^ < 0.10. 



* On behalf of the CDF Collaboration 
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II. DIFFRACTIVE DUET PRODUCTION 



The gluon and quark content of the interacting partons can be investigated by comparing 
single diffractive (SD) and non diffractive (ND) events. SD events are triggered on a leading 
anti-proton in the RPS and at least one jet, while the ND trigger requires at least one jet in the 
calorimeters. The ratio of SD to ND dijet production rates (Njj) is proportional to the ratio of 
the corresponding structure functions (Fjj), RsD_(x,f,t) = ^/ '^^2^' ~ " ^, , and can 
be measured as a function of the Bjorken scaling variable x = XBj[4]. In the ratio, jet energy 
corrections approximately cancel out, thus avoiding dependence on Monte Carlo (MC) simulation. 
Diffractive dijet rates are suppressed by a factor of 0(10) with respect to expectations based on 
the proton PDF obtained from diffractive deep inelastic scattering at the HERA ep collider 
The SD/ND ratios (i.e. gap fractions) of dijets, W, b-quark, J/ip production are all approximately 
1%, indicating that the suppression factor is the same for all processes and it is related to the gap 
formation. 

In Run II, the jet Et spectrum extends to Ej^^ ^ 100 GeV, and results are consistent with those 
of Run I[5], hence confirming a breakdown of factorization. Preliminary results indicate that the 
ratio does not strongly depend on E^ = Q'^ in the range 100 < < 10,000 GeV^ (Fig. [U left). 
The relative normalization uncertainty cancels out in the ratio, and the results indicate that the 
evolution, mostly sensitive to the gluon density, is similar for the proton and the Pomeron. A novel 
technique [g] to align the RPS is used to measure the diffractive dijet cross section as a function of 
the t-slope in the range up to ~ 4, 500 GeV/c^ (Fig. [U right). The shape of the t distribution 
does not depend on the value, in the region < |t| < 1 GeV^. Moreover, the \t\ distributions 
do not show diffractive minima, which could be caused by the interference of imaginary and real 
parts of the interacting partons. 



III. DIFFRACTIVE W/Z PRODUCTION 

Studies of diffractive production of the W/Z bosons are an additional handle to the understand- 
ing of diffractive interactions. At leading order (LO) diffractive W/Z bosons are produced by a 
quark interaction in the Pomeron. Production through a gluon can take place at NLO, which is 
suppressed by a factor and can be distinguished by the presence of one additional jet. 

In Run I, the CDF experiment measured a diffractive W boson event rate Rw = 1.15 it 
0.51 (stat)±0.20 (syst)%. Combining the Rw measurement with the dijet production event rate 
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FIG. 1: Left: Ratio of diffractive to non-diffractive dijet event rates as a function of xsj (momentum fraction 
of struck parton in the anti-proton) for different values of Et'^ = Q^; Right: Measured |t|-distributions for 
soft and hard diffractive events. 



(which takes place both through quarks and gluons) and with the b-production rate allows the de- 
termination of the gluon fraction carried by the Pomeron which can be estimated to be 54^1^% Q]- 
In Run II, the RPS provides an accurate measurement of the fractional energy loss (^) of the 
leading hadron (Fig. [21 left), removing the ambiguity of the gap survival probability. The innovative 
approach of the analysis takes advantage of the full W ^ Iv event kinematics including the neutrino. 
The missing transverse energy (^t) is calculated as usual from all calorimeter towers, and the neu- 
trino direction (i.e. r/j^) is obtained from the comparison between the fractional energy loss measured 
in the Roman Pot spectrometer (C^^^) and the same value estimated from the calorimeters (^'^"'): 

^RPS_^cai ^Sr.e-^-. The reconstructed W mass (Fig. El right) yields Mw = 80.9 ±0.7 GeV/c^, 

J 

in good agreement with the world average value oiMw = 80.398±0.025 GeV/c^[9]. After applying 
the corrections due to the RPS acceptance, trigger and track reconstruction efficiencies, and taking 
into account the effect of multiple interactions, both W and Z diffractive event rates are calculated: 
Rw = 0.97±0.05(stat) ± 0.10(syst)%, and i?z=0.85±0.20 (stat)±0.08 (syst)% |8|. 

Search for exclusive Z candidate events, i.e. pp ^ p + Z + p, has been explored with a null 
result. In the SM, the process takes places through photo-production. The search requires that 
nothing else is found in the detector, except the two leptons from Z — )■ //. The method consists in 
comparing the total energy in the calorimeter (Mx) to the dilepton invariant mass (Mu). Exclusive 
events are expected to be found on the diagonal Mx=Mii. Some effects that may artificially change 
the "exclusive" behavior. By increasing the calorimeter thresholds the value of Mx moves closer to 
the diagonal Mx=Ma. Because of charge conservation, W bosons cannot be produced exclusively 
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FIG. 2: Calorimeter i^'^"' distribution in W events with a reconstructed Roman Pot track {left). Due to 
the neutrino, ^^'^^ < ^^^^ is expected. The difference S,^^^ — £,''°'^ is used to determine the W boson mass 
{right). 



and are used as the control sample. Additional control of the background is performed by looking 
at "empty crossings", where no tracks are reconstructed and calorimeter noise is the dominant 
effect. No exclusive candidate are found in the data. 



IV. FORWARD JETS 



An interesting process is dijet production in double diffractive (DD) dissociation. DD events 
are characterized by the presence of a large central rapidity gap and are presumed to be due to 
the exchange of a color singlet state with vacuum quantum numbers. A study of the dependence 
of the event rate on the width of the gap was performed using Run I data with small statistics. 
In Run II larger samples are available. Typical luminosities (£ ~ 1 10 x lO'^^cm^^sec"^) during 
normal Run II run conditions hamper the study of gap "formation" due to multiple interactions 
which effectively "kill" the gap signature. Central rapidity gap production was studied in soft 
and hard diffractive events collected during a special low luminosity run {C ~ 10^^ cm~^sec~^). 
Figure [3] (left) shows a comparison of the gap fraction rates, as function of the gap width (i.e. A77) 
for minimum bias (MinBias), and MP jet events. Event rate fraction is calculated as the ratio 
of the number of events in a given rapidity gap region divided by all events: Rgap = Ngap/Naii- 
The fraction is approximately 10% in soft diffractive events, and approximately 1% in jet events. 
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Shapes are similar for both soft and hard processes, and gap fraction rates decrease with increasing 
Arj. The MP jets of gap events are produced back-to-back (Fig. [3l right). 
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FIG. 3; Left: Event rate gap fraction defined as Rgap = Ngap/Naii, for minimum bias (MinBias) and MP 
jet events with Et > 2(4) GeV; Right: Azimuthal angle difference distribution of the two leading jets 
in a DD event with a central rapidity gap (|?7^°^| < 3.3). 



V. EXCLUSIVE PRODUCTION 

The first observation of the process of exclusive dijet production can be used as a benchmark 
to establish pred ictions on exclusive diffractive Higgs production, a process with a much smaller 



cross section 
exclusive di 
production 



lot]. A wide range of predictions was attempted to estimate the cross section for 



et and Higgs production. In Run I, the CDF experiment set a limit on exclusive jet 



111 ]. First observation of this process was made in Run II. The search strategy is based 
on measuring the dijet mass fraction (Rjj), defined as the ratio of the two leading jet invariant mass 
divided by the total mass calculated using all calorimeter towers. An exclusive signal is expected 
to appear at large Rjj values (Fig. [H left). The method used to extract the exclusive signal from 
the Rjj distribution is based on fitting the data to MC simulations. The quark/gluon composition 
of dijet final states can be exploited to provide additional hints on exclusive dijet production. The 
Rjj distribution can be constructed using inclusive or b-tagged dijet events. In the latter case, as 
the gg — t- qq is strongly suppressed for rriq/M'^ — t- {Jz = selection rule), only gluon jets will be 
produced exclusively and heavy flavor jet production is suppressed. Figure H] (center) illustrates 
the method that was used to determine the heavy-fiavor composition of the final sample. The 
falling distribution at large values of Rjj (Rjj > 0.7) indicates the suppression of the exclusive 
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events. The CDF result favors the model in Ref. [1^ (Fig. HI right). Details can be found in 
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FIG. 4: Left: Dijet mass fraction Rjj in inclusive DPE dijet data. An excess over predictions at large Rjj 
is observed as a signal of exclusive dijet production; Center: Ratio of b-tagged jets to all inclusive jets as a 
function of the mass fraction Rjj. The error band corresponds to the overall systematic uncertainty; Right: 
The cross section for events with Rjj > 0.8 is compared to predictions. 



Exclusive e~^e~ and di-photon production were studied using a trigger that requires forward 
gaps on both sides of the interaction point and at least two energy clusters in the electromagnetic 
calorimeters with transverse energy Et > 5 GeV. All other calorimeter towers are required to be 
below threshold. In the di-electron event selection, the two tracks pointing at the energy clusters are 



li- A 



allowed. The CDF experiment reported the first observation of exclusive e~^e~ production 
total of 16 77 — )• e~^e~ candidate events are observed, consistent with QED expectations. Exclusive 
di-photon events can be produced through the process gg — ?• 77. Three candidate events were 
selected, where one is expected from background sources (i.e. vr'^vr*^). A 95%C.L. cross section limit 
of 410 pb can be set |l5|, about ten times larger than expectations [l^ . 



VI. CONCLUSIONS 



The results obtained during the past two decades have led the way to the identification of striking 
characteristics in diffraction. Moreover, they have significantly contributed to an understanding 
of diffraction in terms of the underlying inclusive parton distribution functions. The regularities 
found in the Tevatron data and the interpretations of the measurements can be extrapolated to 
the LHC era. At the LHC, the diffractive Higgs can be studied but not without challenges, as 
triggering and event acceptance will be difficult. Still, future research at the Tevatron and at the 
LHC holds much promise for further understanding of diffractive processes. 
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Forward Physics at LHCb 
Prospects for the Study of Diffractive Interactions- 
Michael Schmelling^' 

^ MPI for Nuclear Physics, Saupfercheckweg 1, D-691 17 Heidelberg, German^ 

LHCb, the smaUest of the large LHC experiments, is a forward spectrometer covering 
the angular range 2 < r/ < 5 with tracking, calorimetry and particle identihcation. Partial 
coverage of the backward hemisphere is also provided by the vertex detector (Vertex Locator, 
VeLo), a silicon strip detector surrounding the interaction region. Generator level Monte 
Carlo studies suggest that using the VeLo to ask for a rapidity gap of A77 — 2.5 in the 
backward region allows to select event samples dominated by diffractive processes. Making 
use of the excellent tracking, vertexing and particle identihcation capabilities of the LHCb 
detector, the characteristics of diffractive particle production thus can be studied in detail 
in the forward acceptance covered by the experiment. 







I. INTRODUCTION 



With the startup of LHC an energy regime has become experimentally accessible which will 
allow to probe fundamental physics with unprecedented sensitivity. Although built with the focus 
on finding the Higgs boson and doing searches for physics beyond the standard model, also basic 
questions of particle production in high energy collisions have to be addressed. Minimum bias 
physics at the LHC differs from previous studies at hadron machines in that the center-of-mass 
energy has reached a level, where even collisions between very soft partons can contribute to final 
state particle production. The typical scale is given by the requirement xi ■ X2 • s > 4m^, i.e. 
X ~ 2m^/-y/s. Since the parton densities at small x and small momentum transfer rise faster than 
1/x, it is expected that multi-parton interactions become important and that a new holistic picture 
including diffractive processes for describing such interactions is required. In this paper some 
thoughts addressing this issue are discussed together with first ideas how the LHCb experiment 
can contribute in this area. 



*on behalf of the LHCb Collaboration 
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II. THE LHCB EXPERIMENT 

The LHCb detector [l| is constructed as a forward spectrometer, covering the angular range of 
15 < ^ < 300mrad with respect to the beam axis. A schematic view of the experiment is shown 
in Fig.[TJ Momentum measurement is performed with a dipole magnet with a field integral of 
4 Tm. In front of the magnet the Vertex Locator (VeLo) surrounds the interaction region. Going 
downstream, a first RICH detector and the so-called TT tracking station are still located in front 
of the magnet. Immediately behind the magnet follows the second part of the tracking system, 
consisting of a high granularity Inner Tracker (IT) in the region of large particle densities close to 
the beam pipe and the Outer Tracker system at larger transverse distances. VeLo, TT and IT are 
silicon strip detectors, the OT consists of straw tubes. Following the tracking system is a second 
RICH detector, a pre-shower and scintillating pad detector (SPD/PS), electromagnetic calorime- 
ter (ECAL), hadron calorimeter (HCAL) and muon system for the identification of electrons and 
photons, neutral hadrons and muons, respectively. The RICH detectors allow pion, kaon, proton 
separation in the momentum range between 2 <p < lOOGeV/c. The detector is constructed such 
that it offers tracking, calorimetry and particle identification over most of its forward acceptance. 




FIG. 1: Schematic view of the LHCb single arm forward spectrometer. The interaction region is located on 
the left inside the Vertex Locator. The tracking system and the RICH detectors for particle identification are 
installed both before and after the dipole magnet, calorimetry and the muon system are located downstream 
of the magnet. 

The Vertex Locator has 21 double-layer sensor planes around the interaction region for measur- 
ing space points, plus two additional planes providing only radial track coordinates. The layout of 



100 



the VeLo is shown in Fig.[2l It has a larger angular acceptance than the rest of the tracking system 
and covers even part of the backward hemisphere. However, being located outside of the magnetic 
field, VeLo track segments do not have momentum information. Furthermore, since at least three 
planes are required to reconstruct a track segment, the VeLo is blind in the central region. Charged 
particle tracks are reconstructed in the rapidity ranges —4 < r] < —1.5 and 1.5 < rj < 5. As will be 
shown below, the large angular coverage of the VeLo is vital for the study of diffractive processes. 



1 m 




mrad 



Interaction region o = 5.3cm 

FIG. 2: Layout of the LHCb Vertex Locator (VeLo) in the horizontal (x, z)-plane, with the z-axis along the 
direction of the proton beams. 21 sensor planes measure space points, the two most backward (— z) layers 
provide only radial coordinates of charged particle tracks. 



III. INELASTIC PROTON-PROTON INTERACTIONS 

A schematic view of different types of inelastic pp- interactions is presented in Fig.[3l Here the 
basic distinction is colour(-octet)- and colour-singlet exchange, respectively, between the colliding 
protons. Colour-exchange implies that the structure of both protons is resolved with the conse- 
quence that the colour fields stretched between the partons lead to particle production in the full 
rapidity range. In contrast, colour-singlet exchange is phenomenologically described by pomerons 
coupling to the protons as a whole. No colour is transferred and the protons can either scatter 
elastically or be excited into a high mass state which then decays to produce multi-particle final 
states. Depending on whether only one or both protons are excited these processes are referred to as 
single- or double-diffractive scattering. An example for a higher order process involving pomerons 
is double pomeron exchange, where both protons stay intact and the two pomerons interact to 
form a massive central system. 

The above classification of pp-inter actions is most adequate at small momentum transfers. In the 
language of QCD then colour exchange means gluon exchange and the pomeron can be understood 
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FIG. 3: Schematic classification of inelastic proton-proton collisions. The upper row shows some born-level 
type diagrams for different classes of interactions, the lower row illustrates the angular range into which 
particles produced in the collision are emitted. Note that in the upper row rapidity runs from top to bottom, 
while it goes from left to right in the lower row. Also, while the born-level type diagrams represent scattering 
amplitudes, the particle flow corresponds to cross sections, i.e. the squares of the respective amplitudes. 



as a colour-singlet two gluon state. With increasing momentum transfer, however, the simple 
picture breaks down. The two protons are resolved into an increasing number of partons and 
the interaction is described by ladder-diagrams of all possible topologies. Diffractive and non- 
diffractive scattering is no longer an unambiguous classification, and even the notion of colour 
singlet exchange becomes frame dependent ^. It follows that eventually a unified description of 
hadron-hadron collisions is required which covers diffractive and non-diffractive physics within a 
coherent framework Q]. 

FigureO illustrates another conceptual problem with respect to diffractive and non-diffractive 
interactions. While the upper row represents Born-level type diagrams of the different processes, 
i.e. amplitudes contributing to the inelastic proton-proton scattering, the lower one is a pictorial 
representation of the cross-section. In many Monte Carlo models, such as e.g. PYTHIA J|, the 
different components making up the total cross-section are generated independently. Interference 
terms between the different amplitudes are ignored. Physics-wise, however, these terms certainly 
contribute, as it is experimentally not possible to distinguish the fast forward proton from a single- 
diffractive scattering from a non-diffractive interaction where a fast proton is generated in the 
fragmentation process [5]. An unambiguous classification of an event into one of several types of 
inelastic processes thus is impossible. Furthermore, any such separation within the context of a 
specific Monte Carlo simulation is to some extent arbitrary since parameter tuning generally allows 
to trade e.g. a larger diffractive cross-sections against a smaller non-diffractive part by choosing a 
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different setting of the hadronization parameters. 

In the past experimental measurements often were corrected for e.g. single diffractive contribu- 
tions to the cross section. From the previous discussion it is clear, that such attempts to focus on 
the non-diffractive cross section are always model dependent, even though they work to a certain 
degree since, in the language of Fig.[3l some amplitudes are dominant in specific regions of phase 
space. Nevertheless, a better approach would be to avoid any such artificial classifications, and 
instead perform measurements subject to experimental cuts which enhance or suppress certain 
contributions to the cross section. That way the results of the measurements do not rely on a 
particular model, and measurements and their interpretation are cleanly separated. 

IV. MONTE CARLO STUDIES 

To study the prospects for experimentally studying the properties of events with dominantly 
diffractive contributions, a simple generator level study has been performed. The study is based 
on PYTHTIA 8. 135 available from [6]. Single proton-proton collisions with a center-of-mass energy 
^/s = 7TeV were generated with process selection pythia. readStringC'Sof tQCD : all=on") . 

The study focuses on a measurement of the inclusive charged particle transverse momentum 
spectrum for all events and for events with an enhanced fraction of diffractive contributions. For 
this tracking based study the VeLo was simulated with its nominal geometry. A track was assumed 
to be measured by the VeLo if at least three stations were hit. The event selection was based on 
the VeLo track segments only. For accepted events the transverse momentum spectrum then was 
determined using all tracks with a VeLo-segment and within the acceptance of the tracking system 
behind the magnet. The latter was approximated by the requirement p > 2GeV/c the pseudo- 
rapidity range 2 < < 5. 

A diffraction-enriched event sample was selected by exploiting the fact that diffractive events 
are characterized by enhanced probabilities for large rapidity gaps in the final state particle distri- 
bution. To obtain a quantitative measure for the level of enrichment which can be achieved, the 
PYTHIA process type was analyzed for all events. While evidently giving model dependent estimates 
for the fractions of different events, the qualitative picture is expected to be generic. 

Denoting by and np the number of VeLo track segments in the backward (?? < 0) and 
forward (r/ > 0) hemispheres, two selection criteria were studied: (a) ns + np > and (b) 
ub == && np > 0. Criterion (a) corresponds to the so called micro-bias trigger of LHCb, which 
is close to 100% efficient for non-diffractive pp-collisions. The second criterion asks for no charged 
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FIG. 4: Generator level MC study: Mix of single diffractive, double diffractive and non-diffractive events 
generated by PYTHIA 8 . 135. Here SDl refers to single diffractive scattering where the excited proton travels 
in the direction of the LHCb detector, in SD2 the decaying heavy mass moves into the opposite direction. 
The left hand plot shows the mix of events generated (black), passing the micro-bias trigger (red) and the 
diffraction selection (blue). The right hand plot shows the selection efficiencies for the two cases. 



tracks in the backwards acceptance of the VeLo, i.e. it corresponds to a requirement of a rapidity 
gap Ar] = 2.5 for charged tracks. 

Results are shown in Figs.|l]and[5j The left hand plot of Fig.|3]shows the mix of single diffractive, 
double diffractive and non-diffractive events generated by PYTHIA 8. 135. One clearly sees that the 
requirement of a rapidity gap in the backwards region almost completely suppresses non-diffractive 
events while keeping between 20% and 30% of diffractive interactions. The comparison of the 
transverse momentum spectra in Fig. [5] shows very good agreement between generated and observed 
distribution, i.e. a robust measurement comparing the fully inclusive transverse momentum spectra 
and the spectra in events dominated by diffraction seems feasible. Other observables like charge 
ratios, the production cross-sections for identified particles or particle ratios are a natural extension 
of these studies. 



V. SUMMARY AND OUTLOOK 



A model independent approach toward the study of diffractive particle production in minimum 
bias events has been presented. Generator level Monte Carlo studies suggest that asking for a 
rapidity gap of Ar] = 2.5 in the backward region of LHCb VeLO allows to select event samples 
dominated by diffractive processes. Making use of the excellent tracking, vertexing and particle 
identification capabilities of the LHCb detector, the characteristics of particle production in those 
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FIG. 5: Generator level MC study: comparison of generated and observable transverse momentum spectra 
of charged particles in the pseudo rapidity range 2 < r/ < 5 when asking for a micro-bias trigger (left) or 
a rapidity gap (right). The spectra are normalized to the number of accepted events. Finite resolution or 
imperfect tracking efhciency has not been modeled. For the diffraction-enhanced event sample the MC-truth 
level has been defined to include only events of type SDl and DD. The losses at small transverse momenta 
are due to incomplete geometric coverage at low px and low 77. 



events can be studied in detail in the pseudo-rapidity range 2 < r/ < 5. 
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The forward physics program of the CMS experiment at the LHC spans a broad range 
of diverse physics topics including studies of low-x QCD and diffractive scattering, muhi- 
parton interactions and underlying event structure, 7-mediated processes and luminosity 
determination, Monte Carlo tuning and even MSSM Higgs discovery in central exclusive 
production. In this article, the forward detector instrumentation around the CMS interaction 
point is described and the prospects for diffractive and forward physics using the CMS 
forward detectors are summarized. In addition, first observation of forward jets as well as 
early measurements of the forward energy flow in the pseudorapidity range 3.15 < \ri\ < 4.9 
at = 0.9 TeV, 2.36 TeV and 7 TeV are presented. 
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I. THE CMS EXPERIMENT AT THE LHC 

The Compact Muon Solenoid (CMS) [1] is one of two general-purpose particle physics detectors 
built at the Large Hadron Collider (LHC) at CERN. The detector has been designed to study 
various aspects of proton-proton (pp) collisions at ^/s =14 TeV and heavy-ion(Pb-Pb) collisions 
at ^/s =5.5 TeV, that will be provided by the LHC at a design luminosity of lO'^'* cm~^s~^ and 
of 10^'' cm~^s~^, correspondingly. To enhance the physics reach of the experiment the CMS 
subcomponents must provide high-precision measurements of the momentum and the energy of 
collision-products. The CMS detector comprises the tracking system covering the pseudorapidyty 
range —2.5 < rj < 2.5 and the calorimetry system covering the pseudorapity range — 5 < r/ < 5. 
In addition to that, CMS includes several very forward calorimeters, whose design and physics 
potential will be described later in this article. It should be emphasized that the CMS detector 
is one of the largest scientific instruments ever built. It comprises about 76.5 millions of readout 
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channels in total. The detector has been designed, constructed and currently operated by the 
collaboration consisting of more than 3500 scientists from 38 countries. 

First collision data taking at CMS took place in November 2009. Since then and by the end of 
May 2010, CMS has collected around 10 nb~^ of collision data. It should be noted that the quality 
of collected data is rather good: more than 99% of CMS readout channels are operational and the 
CMS data taking efficiency is above 90%. Several tens of pb~^ of the pp collision data are expected 
to be collected by the end of 2010. 

II. FORWARD DETECTORS AROUND THE CMS INTERACTION POINT 

The maximum possible rapidity at the LHC in pp collisions at s/s = 14 TeV is Umax = 
In ("v/i/mTr) ~ 11.5 and one of the great features of the CMS experiment is that it includes 
several subdetectors covering the kinematic region at very small polar angles and so, large values 
of rapidity. A schematic view of CMS forward detectors is shown in Figure 1. As can be seen, the 
CMS forward instrumentation consists of the Hadronic Forward calorimeter (HF), the CASTOR 
and ZDC calorimeters. All of them are sampling calorimeters. That is, they are made of repeating 
layers of a dense absorber and tiles of scintillator. A separate experiment TOTEM as well as proton 
detectors FP420 are additional forward detectors around the CMS interaction point (IP5). They 
further extend the forward reach available around IP5. 

A. HF 

The CMS HF detector [2] includes two calorimeters HF+ and HF-, which are located at a 
distance of 11.2 m on the both sides from the IP5 covering the pseudorapidity range 3 < |r7| < 5. 
The detector is designed to carry out the measurements of the forward energy flow and forward jets. 
The HF is a Cerenkov sampling calorimeter which uses radiation hard quartz fibers as the active 
material and steel plates as the absorber. The signal in the HF is produced when charged shower 
particles pass through the quartz fibers with the energy above the Cerenkov threshold (190 keV 
for electrons). The generated Cerenkov light is then collected by air-code light guides, which are 
connected to photo-multipliers tubes PMTs. The detector fibers run parallel to the beamline and 
are bundled to form 0.175 x 1.175 (A?7 x A^) towers. Half of the fibers run over the full depth of the 
absorber, whereas the other half starts at a depth of 22 cm from the front of the detector. These 
two sets of fibers are read out separately. Such a structure allows to distinguish showers generated 
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FIG. 1: Layout of the forward detectors around the CMS interaction point. 



by electrons and photons, which deposit a large fraction of their energy in the first 22 cm, from 
those generated by hadrons, which produce signals in both calorimeter segments. The detector is 
embedded into a cylindrical steel structure with the outer radius of 131 cm and the inner radius 
of 12.5 cm to accommodate the beam pipe. It is azimuthally subdivided into 20" modular wedges, 
each of which consists of two azimuthal sectors of 10*^. The detector extends over 10 interaction 
lengths and includes 1200 towers in total. 



B. CASTOR 



The CASTOR (Centauro And STrange Object Reseacrh) detector [3] is a quartz-tungsten 
Cerenkov sampling calorimeter, which is located at a distance of 14.4 m from the IP5 and covering 
the pseudorapidity range —6.6 < r] < —5.2. The tungsten metal has been chosen as the absorber 
medium in CASTOR, since it provides the smallest possible shower size. In this detector, the 
radiation hard quartz plates used as the active medium are tilted at 45*^ to efficiently capture the 
Cerenkov light produced by relativistic particles passing the detector. As in the case of the HF, the 
produced Cerenkov light is collected by air-code light guides that are connected to PMTs, which 
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produce signals proportional to the amount of light collected. The CASTOR detector is a compact 
calorimeter with the physical size of about 65 cm x 36 cm x 150 cm and having no segmentation in 
rj. It is embedded into a skeleton, which is made of stainless steel. The detector consists of 14 lon- 
gitudinal modules, each of which comprises 16 azimuthal sectors that are mechanically organized 
in two half calorimeters. First 2 longitudinal modules form the electromagnetic section, while the 
other 12 modules form the hadronic section. In the electromagnetic section, the thicknesses of the 
tungsten and quartz plates are 5.0 and 2.0 mm respectively, whereas in the hadronic section the 
corresponding thicknesses are 10.0 and 4.0 mm. With this design, the diameter of the showers 
of electrons and positrons produced by hadrons is about one cm, which is an order of magnitude 
smaller than in other types of calorimeters. The detector has a total depth of 10.3 interaction 
lengths and includes 224 readout channels. 



C. ZDC 



The CMS ZDC (Zero Degree Calorimeter) detector [4] consists of two calorimeters that are 
located inside the TAN absorbers at the ends of the straight section of the LHC beampipe at a 
distance of 140 m on both sides from the IP5. These are Cerenkov sampling calorimeters with 
quarz fibers as the active material and tungsten plates as the absorber material. The ZDC detec- 
tor is designed to measure neutrons and very forward photons providing detection coverage in the 
pseudorapidity region \r]\ > 8.4. Each ZDC is made up of separate electromagnetic and hadronic 
sections. The electromagnetic section consists of 33 layers of 2 mm thick tungsten plates and 33 
layers of 0.7 mm thick quartz fibers. The hadronic section is made of 24 layers of 15.5 mm thick 
tungsten plates and 24 layers of 0.7 mm thick quartz fibers. The electromagnetic section is seg- 
mented into 5 horizontal individual readout towers, whereas the hadronic section is longitudinally 
segmented into 4 readout segments. The tungsten plates are oriented vertically in the electromag- 
netic section whereas they are tilted by 45° in the hadronic section. The detector is read out via 
aircore light guides and PMTs. It has a total depth of 6.5 interaction length. 



D. TOTEM and FP420 



TOTEM [5] is an independent experiment at the CMS interaction point whose main objectives 
are the precise measurement of the total pp cross-section and a study of elastic and diffractive scat- 
tering at the LHC. To achieve optimum forward coverage for charged particles, TOTEM comprises 
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two tracking telescopes, Tl and T2, that are installed on both sides from the IPS in the pseudo- 
rapidity region 3.1 < \r]\ < 6.5, and Roman Pot stations that are located at distances of ±147 m 
and ±220 m from the IP5. The Tl telescope is located in front of HF and consists of 5 planes of 
cathode strip chambers, while the T2 telescope is located in front of CASTOR and comprises 10 
planes of gas electron multipliers. For efficient reconstruction of very forward protons, silicon strip 
detectors are housed in the Roman Pot stations. 

FP420 [6] is a proposed detector system, which is supposed to provide proton detection at a 
distance of ±420 m from the IPS. The FP420 detector comprises a silicon tracking system that can 
be moved transversely and measure the spatial position of protons, which have been bent out by the 
LHC magnets due to the loss of a small fraction of their initial momentum. The potential physics 
topics that can be studied with this detector system include Higgs central exclusive production as 
well as a rich QCD and electroweak program. 

III. PHYSICS PROGRAM 

Extending the physics reach of CMS, the program for forward physics includes studies of low- 
X QCD and diffractive scattering, multi-parton interactions and underlying event structure, 7- 
mediated processes and luminosity determination. It is also supposed to contribute to the discovery 
physics via searches of MSSM Higgs in central exclusive production. 

A. Low-x QCD 

A study of QCD processes at a very low parton momentum fraction x = Pparton/Phadron is a 
key to understand the structure of the proton, whose gluon density is poorly known at very low 
values of x. Low-x QCD dynamics can be studied in pp collisions if the parton momentum fraction 
of one of the colliding protons xi is significantly larger than the parton momentum fraction of the 
other colliding proton X2 {xi » ,T2). The result of such a collision is a creation of either jets, 
prompt-7 or Drell-Yan electron pairs at very low polar angles in the very forward region of the 
detector. Low-x QCD studies at CMS will be a continuation of studies of deep inelastic scattering 
in electron-proton collisions at HERA, where low-x QCD dynamics has been explored down to 
values of 10~^. Measurements at HERA have shown that the gluon density in the proton rises 
rapidly with decreasing values of x. As long as the densities are not too high this rise can either be 
described by the DGLAP model [7] that assumes strong ordering in the transverse momentum fer 



110 




log (x ) 

^10^ 1,2' 

FIG. 2: log{xi^2) distribution of two partons producing at least one jet with Et > 20 GeV in the HF 
acceptance. 

or by the BFKL model [8] that assumes strong ordering in x and random walk in kT- Eventually 
at low enough x, the gluon-gluon fusion effects become important saturating the growth of the 
parton densities. 

At the LHC the minimum accessible x in pp collisions decreases by a factor of about 10 for 
each 2 units of rapidity. This implies that a process with a hard scale of Q ~ 10 GeV and within 
the CAST0R/T2 detector acceptance can probe quark densities down x ~ 10^®. Such processes 
include the production of forward jets and Drell-Yan electron pairs. 

1. Forward Jets 

A low-x parton distribution function (PDF) of the proton can be constrained by measuring 
single inclusive jet cross-section in HF. Figure 2 illustrates the log{xi^2) distribution for parton- 
parton scattering in pp collisions at ^/s = 14 TeV requiring at least one jet with the transverse 
energy above 20 GeV in the HF acceptance. As can be seen, by measuring forward jets in HF 
one can probe x values as low as 10~^. A detailed analysis of fully simulated and reconstructed 
QCD jet events generated with PYTHIA in the range px = 20 GeV/c~200 GeV/c in pp collisions 
at ^/s = 14 TeV for an integrated luminosity of 1 pb~^ shows that the momentum resolution 
for forward jets in HF is about 18% at px = 20 GeV/c and is gradually decreasing to 12% at 
PT > 100 GeV/c [9]. 
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A possibility to gain information on the full QCD evolution to study high order QCD reactions 
can be provided by measuring forward jets in the CASTOR calorimeter, that will allow to probe 
the parton densities down 10~^. Apart from that, it has been found that a BFKL like simulation 
predicts more hard jets in the CASTOR acceptance than the DGLAP model. Therefore, measure- 
ments of forward jets in CASTOR can be used as a good tool to distinguish between DGLAP and 
non-DGLAP type of QCD evolution. 

Further studies of low- a; QCD can be made with Mueller-Navalet dijet events, which are char- 
acterized by two jets with similar pT but large rapidity separation. By measuring Mueller-Navalet 
dijets in CASTOR and HF one can probe BFKL-like dynamics and small-x evolution. 

2. Drell-Yan 

Low-x proton PDFs can also be constructed by measuring electron pairs produced via the 
Drell-Yan process qq ^ j* ^ e+e" within the acceptance of CASTOR and T0TEM-T2 station, 
whose usage is essential for detecting these events. Figure 3 illustrates the distribvition of the 
invariant mass M of the ee system against the parton momentum fraction X2 of one of the quarks, 
where X2 is chosen such that xi » X2- In this figure, the solid line indicates the kinematic limit, 
whereas the region between the dotted lines is the acceptance window for both electrons to be 
detectable in CAST0R/T2. The green points show the events with at least one electron lying 
in CAST0R/T2 acceptance and the blue points indicate the events with both electrons present 
within the CAST0R/T2 acceptance, while the black points correspond to any of the Drell-Yan 
events generated with PYTHIA. As can be seen, by measuring two electrons in the CAST0R/T2 
acceptance one can access x values down to 10~^ for M > 10 GeV [10]. Futhermore, measurements 
of Drell-Yan events in the CAST0R/T2 acceptance can be used to study QCD saturation effects. 
It has been found that the Drell-Yan production cross section is suppressed roughly by a factor of 
2 when using a PDF with saturation effects compared to one without. 

B. Diffraction 

In pp collisions a diffractive process is a reaction pp — > XY, where X and Y can either be 
protons or low-mass systems which may be a resonance or a continuum state. In all cases, the 
final states X and Y acquire the energy approximately equal to that of the incoming protons and 
carry the quantum numbers of the proton as well as are separated by a Large Rapidity Gap (LRG). 
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FIG. 3: Acceptance of the CAST0R/T2 detectors for Drell-Yan electrons. See text for details. 

Diffraction in tlie presence of a liard scale can be described with perturbative QCD by the exchange 
of a colourless state of quarks or gluons, whereas soft diffraction at high energies is described in the 
Regge Theory [11] as a colourless exchange mediated by the Pomeron having the quantum numbers 
of the vacuum. The cross section of hard diffractive processes can be factorized into generalized 
parton distributions and diffractive parton distributions functions (dPDF), which contain a valuable 
information about low-x partons. However, the factorization becomes broken when scattering 
between spectator partons takes place. This effect is quantified by the so-called rapidity gap 
survival probability that can be probed by measuring the ratio of diffractive to inclusive processes 
with the same hard scale. At the Tevatron, the ratio is found to be 0(1%), whereas theoretical 
expectations at the LHC vary from a fraction of a percent to up to 30% [12]. 

The two main types of diffractive processes occurring in pp collisions are the single diffractive 
dissociation (SD) where one of the protons dissociates and the double diffractive dissociation (DD) 
where both protons are scattered into a low-mass system. The single-diffractive productions of W 
and dijets are in particular very interesting processes to study, since they are sensitive to the quark 
and gluon content of the PDFs, correspondingly. They both are hard diffractive processes that 
can provide information on the rapidity gap survival probability. A selection of such events can be 
performed using the multiplicity distributions of tracks in the central tracker and calorimeter towers 
in HF plus CASTOR exploiting the fact that diffractive events on average have lower multiplicity 
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in the central region and in the "gap side" than non-difractive ones. Feasibihty studies to detect 
the SD productions of W [13] and dijets [14] have shown that the diffractive events peak in the 
regions of no activity in HF and CASTOR. 

C. Exclusive dilepton production 

Another interesting topic that is going to be studied at CMS is the exclusive dilepton production 
pp —7- ppl^l^ , which can either occur via T photoproduction 7p — > T — ?► l^l^ or via the pure QED 
process 77 — t- l^l^ that has been observed by the CDF experiment at the Tcvatron [15]. The 
latter is an elastic process whose production cross section is precisely known. As a result, it 
can potentially serve as an ideal calibration channel and is going to be used for measuring the 
luminosity at the LHC. Using this process an absolute luminosity calibration with the accuracy 
of 4% is feasible with 100 pb~^ of data [16]. The dominant background source for this mode is 
inelastic processes, where one of the proton in the process does not stay intact but dissociates. It 
can be significantly suppressed by applying a veto condition on activity in CASTOR and ZDC. 
Exclusive dilepton production occurring via T photoproduction is also a mode of interest, since 
the cross section of the T photoproduction process is sensitive to the generaliszed PDF for gluons 
in the proton. Finally, it should be noted that exclusive dimuon production is an ideal alignment 
channel for the proposed FP420 proton detectors. 

D. Multi-parton interactions and forward energy flow 

Multi-parton interactions (MPI) arise in the region of small-x where parton densities are large 
so that the likelihood of more than one parton interaction per event is high. According to all 
QCD models, the larger the collision energy the greater the contribution from multiple parton 
interactions to the hard scattering process. However, the dependence of the MPI cross section on 
the collision energy is not well known and needs to be studied. A good way to study multiple 
parton interactions is provided by the energy flow in the forward region, which is directly sensitive 
to the amount of parton radiation and MPI. Measurements of the forward energy flow will allow to 
discriminate between different MPI models, which vary quite a lot, and provide additional input to 
the determination of the parameters of the existing MPI models. Furthermore, measurements of 
forward particle production in pp and Pb-Pb collisions at LHC energies should help to significantly 
improve the existing constraints on ultra-high energy cosmic ray models. The primary energy 
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FIG. 4: Display of an event with two forward jets. 

and composition of the ultra-high energy cosmic rays are currently determined from Monte Carlo 
simulations using Regge-Gribov-based approaches [17] (where the primary particle production is 
dominated by forward and soft QCD interactions) with parameters constrained by the existing 
collider data at the Eiab < 10^^ eV, whereas the measured energies of the ultra-high energy cosmic 
rays extend up to 10^*^ eV and even beyond. At the LHC energy of Eiab = 10^'^ eV, a more reliable 
determination of the cosmic ray energy and composition becomes possible. Finally, it should be 
emphasized that the forward energy flow has never previously been measured at a hadron collider. 

IV. FIRST RESULTS FROM CMS 
A. Observation of forward jets 

A search for forward jets in the pseudorapidity range 3 < |?7| < 5 has been made as soon as the 
CMS detector has started to take collision data [18]. One of the first candidates of a forward dijet 
event recorded by CMS at ^/s = 0.9 TeV is shown in Figure 4. The displayed event includes one 
forward jet and one backward jet both with a corrected pT above 10 GeV/c. 
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B. Measurement of the forward energy flow 

Early measurements of the energy flow in the forward region of the CMS detector have been 
made with minimum bias events using the pp collision data sets collected at -y/s = 0.9 TeV and 
2.36 TeV in the fall of 2009 and sX^ = l TeV in March 2010 [19]. To select the events of interest 
the following conditions were imposed. 

First, the Beam Scintillator Counters (BSC) and the Beam Pick-up Timing for the experi- 
ments (BPTX), both are elements of the CMS detector monitoring system, were used to trigger 
the detector readout. The BSC devices are located at a distance of 10.86 m on both sides from 
the interaction point covering the pseudorapidity range 3.23 < < 4.65 and providing hit and 
coincidence signals with a time resolution of about 3 ns. Each BSC comprises 16 scintillator tiles. 
The two BPTX elements are located around the beam pipe at a distance of ±175 m from the 
interaction point providing precise information on the bunch structure and timing of the incoming 
beam with a time resolution better than 0.2 ns. To select the minimum bias events with activity 
in the forward regions, the coincidence between a trigger signal in the BSC scintillators and BPTX 
signals was required for both beams. 

Next, to ensure that the selected event is a collision candidate, the events were required to have 
at least one primary vertex reconstructed from at least 3 tracks with a z distance to the interaction 
point below 15 cm and a transverse distance from the z-axis smaller than 2 cm. Further cuts were 
applied to reject beam-scrapping and beam-halo events. Finally, the energy threshold of 4 GeV 
has been imposed to suppress electronic noise in HF. 

In this study, the measurement of energy flow has been made at detector level in the pseudora- 
pidity range 3.15 < \rj\ < 4.9 covered by the HF calorimeters. The energy flow ratio, estimated in 
this analysis, is defined as 

^Efiow - 1 ' y^) 

dr, 

where is the number of selected events, dE ^ is the energy deposition integrated over (p 

in the region dr], -^/si refers to either 2.36 TeV or 7 TeV, whereas refers to 0.9 TeV. The 
pseudorapidity range is divided into five bins with a size of 0.35 in units of r/ following the transverse 
segmentation of the HF calorimeters. In Figures 5 and 6, the energy flow ratio is shown for different 
collision energies as the average of the HF(-|-) and HF(-) responses. In these plots, uncorrected 
data without systematic uncertainties are compared to simulated events obtained from PYTHIA 
tune D6T. As can be clearly seen, the energy flow gets larger at forward rapidities and with 
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FIG. 5: Energy flow ratio for Vs]" = 2.36 TcV to Vs2 = 0.9 TeV as a function of 77. See text for details. 
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FIG. 6: Energy flow ratio for y/si = 7 TeV to y/s2 = 0.9 TeV as a function of r]. See text for details. 



increasing centre-of-mass energy. Apart from that, it should be noted that the obtained results do 
approximately agree with the Monte Carlo predictions. However, no conclusions on the quality of 
the description can be drawn in this early study due to the missing systematic effects. 
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V. CONCLUSIONS 

A very rich forward physics program can be made with the CMS detector at the LHC due to 
the unprecedented kinematic coverage of the forward region. All the CMS forward detectors have 
been successfully commissioned in 2009 and currently take collision data. The first measurement 
of the forward energy flow has been performed and forward jets at \r)\ > 3 have been observed for 
the first time at hadron colliders. 
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